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ABSTRACT 

Sikeri,  Oddvar,  M.S.C.E.  Purdue  University,  January  1964.  "Hydraulics 
of  Submerged  Arch  Bridges".  Major  Professor:  Dr.  J.  W.  Delleur. 

This  report  is  part  of  a  general  program  sponsored  by  the  State 
Highway  Department  of  Indiana  and  the  U.  S.  Bureau  of  Public  Roads  at 
Purdue  University  on  the  Hydraulics  of  River  Flow  Under  Arch  Bridges. 

The  project  has  been  conducted  since  1958  and  several  reports  on 

different  phases  of  the  study  have  been  submitted.  In  1958  Mr.  S.  T. 

1  * 
Husain   reported  his  work,  "Preliminary  Model  Investigation  of  Hydraulic 

Characteristics  of  River  Flow  Under  Arch  Bridges".  This  was  an  investi- 
gation in  a  12  foot  long,  6  inches  wide  tilting  flume  with  smooth 
boundaries  which  led  to  the  building  of  a  6b   foot  by  5  foot  by  2  foot 

all  steel  tilting  flume.  The  construction  of  the  6b   foot  flume  was 

2 
undertaken  by  Mr.  H.  J.  Owen  and  the  work  was  reported  in  his  masters 

thesis  entitled  "Design  and  Construction  of  Hydraulic  Flume  and  Back- 
water Effects  of  Semi-circular  Constrictions  in < Smooth  Rectangular 
Channel".  Mr.  A.  A.  Sookjr  derived  an  equation  for  the  discharge 
through  a  sharp  edged  semi-circular  constriction  for  the  case  of  free 
surface  flow  through  the  opening.  He  also  extended  Mr.  S.  T.  Husain1 s 
preliminary  studies  to  include  the  effect  of  rough  boundaries.  The 

work  of  Husain,  Sooky  and  Owen  was  summarized  in  a  report  entitled, 

„  3 
"Hydraulics  of  River  Flow  Under  Arch  Bridges  -  A  Progress  Report", 


*  Superscripts  refer  to  references  in  the  bibliography. 
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.  it, 

In  January  1 961  Mr.  P.  F.  Biery  presented  a  report  on  "Hydraulics 

of  Single  Span  Arch  Bridge  Constrictions"  as  part  of  his  work  toward 

a  Master's  degree.     Biery  investigated  the  single  span  arch  bridge 

with  free  surface  flow  through  the  constriction.  The  study  did  not 

include  skew,  eccentricity  or  entrance  rounding.  A  backwater  equation 

was  presented  together  with  a  bridge  design  procedure  for  free  surface 

flow  through  the  opening. 

The  present  report  on  "Hydraulics  of  Submerged  Arch  Bridges" 
covers  the  work  done  by  the  author  from  September  1 962  to  December 
1963.  This  work  was  also  conducted  in  the  large  64  foot  by  5  foot  by 
2  foot  tilting  steel  flume  in  the  Hydraulics  Laboratory  at  Purdue 
University. 

The  existing  instrumentation  was  improved.  The  measurement  of 
flow  velocity  with  the  Prandtl  tube  was  facilitated  by  the  replacement 
of  the  usual  fluid  manometer  by  two  pressure  transducers  and  electric 
accessories.  This  gave  a  more  sensitive  reading  of  small  velocities 
and  a  much  faster  response  to  small  changes  in  the  pressure  differentials, 

A  theoretical  discharge  equation  for  orifice  type  flow  through 
semi-circular  bridge  openings  was  developed.  An  equation  for  the 
head  loss  due  to  the  bridge  constriction  was  also  derived.  Model 
studies  for  submerged  bridge  constrictions  were  made  both  for  smooth 
and  rough  channel  boundaries.  Models  of  two  and  three  dimensions  with 
and  without  skew,  with  one  and  two  openings,  but  having  no  eccentricity 
or  entrance  rounding  were  investigated. 

Through  the  theoretical  and  experimental  analysis  the  following 


XV 


results  were  obtained: 

a)  A  dimensionless  relationship  between  discharge  and 
backwater  elevation. 

b)  A  semi-empirical  expression  to  distinguish  between  the 
occurrence  of  free  surface  flow  and  orifice  type  flow 
through  the  bridge  constriction. 

c)  Two  types  of  graphical  presentation  of  the  test  results 
for  each  geometry.  In  one  presentation,  the  discharge  is 
shown  as  a  function  of  the  ratio  of  the  backwater  depth 
to  the  arch  span,  in  the  other  the  ratio  of  the  backwater 
to  the  normal  depth  is  shown  as  a  function  of  the  Froude 
number  of  the  unconstricted  flow  and  of  the  channel 
opening  ratio. 

d)  Empirical  generalized  backwater  equations  derived  for 
the  several  geometries  from  the  second  presentation 
discussed  under  c. 

e)  Graphs  giving  the  value  of  the  constriction  head  loss 
coefficient  as  a  function  of  the  channel  opening  ratio 
for  different  test  geometries. 

Finally,  suggestions  are  made  for  further  study  of  the  submerged 
arch  bridge  constrictions. 
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.  CHAPTER  I 

Introduction 

How  many  arch  bridges  are  there  in  the  entire  United  States? 
Although  this  figure  is  probably  unknown,  Indiana  alone  has  about 
900  arch  bridges.  By  increasing  the  backwater  upstream  from  the  bridge 
site,  many  of  these  bridge  constrictions  cause  additional  flooding 
during  high  flows;  people  and  valuable  property  are  endangered.  In 
recent  years,  the  problem  of  minimizing  the  backwater  effect  due  to 
bridge  constriction  has  become  increasingly  important.  The  highway 
engineers  are  faced  with  a  multisided  problem;  not  only  do  they  have 
to  build  bridges  to  convey  a  specified  traffic  volume  and  are  safe 
with  respect  to  flooding,  but  they  also  must  find  the  most  economic 
design  possible.  It  is  common  knowledge  to  the  highway  engineer  that 
a  bridge  crossing  will  interfere  with  the  natural  flow  of  the  stream 
and  will  result  in  a  rise  of  stage  upstream  from  the  bridge  and  in 
an  increase  in  velocity  through  the  bridge  constriction.  The  engineers 
are  also  aware  of  alternatives  of  the  inexpensive  soil  fill  compared 
to  a  structural  span.  The  optimum  design,  is  the  shortest  span  that 
will  not  cause  damage  due  to  stage  increase  during  serious  flood 
condition.  In  order  to  meet  these  requirements,  exhaustive  investiga- 
tions, were  undertaken  to  study  the  hydraulic  characteristics  of  the 
water  flow  through  different  bridge  constrictions. 


In  1958,  a  project  was  initiated  in  the  Hydraulics  Laboratory 
at  Purdue  University  in  order  to  study  the  hydraulics  of  stream  flow 
under  arch  bridges.  It  is  sponsored  through  the  Joint  Highway  Research 
Project  in  cooperation  with  U.  S.  Bureau  of  Public  Roads.  Similar 
projects  have  also  been  conducted  at  many  other  universities  through- 
out the  United  States  to  study  different  phases  of  the  hydraulics 
connected  with  flow  through  bridge  crossings. 

The  hydraulics  of  river  flow  through  arch  bridges  is  rather 
complicated,  and  its  investigation  involves  many  variables.  The 
relationships  between  these  variables  cannot  be  handled  completely 
by  means  of  a  theoretical  analysis.  Experimental  models  are  necessary 
in  order  to  investigate  the  problems,  and  even  then  only  the  simplest 
and  most  idealized  conditions  can  be  observed.  In  most  cases  the 
actual  field  conditions  are  so  varied  that  each  one  could  require 
a  separate  model  study.  The  goal  for  this  research  is  to  obtain  some 
empirical  information  about  the  hydraulic  behavior  under  generalized 
conditions,  and  each  field  condition  would  then  have  to  be  given  some 
special  consideration. 

The  present  report  covers  the  investigation  of  the  hydraulics 
of  the  submerged  arch  bridge  constrictions  of  several  geometries. 
The  single-opening,  two-dimensional,  semi-circular  arch  constriction 
of  infinitesimal  length  in  the  flow  direction,  without  skew,  eccentri- 
city or  entrance  rounding  was  studied  first.  Submerged  flow  through 
this  constriction  was  investigated  in  a  smooth  and  in  a  rough  channel 
with  varying  opening  ratios  and  varying  Froude  numbers.  The  stream- 
line pattern,  the  velocity  distribution,  and  the  resulting  rise  in 


the  backwater  upstream  of  the  bridge  constriction,  were  studied. 
The  coefficient  of  contraction,  the  coefficient  of  discharge,  and 
the  coefficient  of  velocity,  were  also  investigated.  The  single 
span  three-dimensional  arch  constriction  (without  skew,  eccentricity 
or  rounding)  of  several  length  simultates  the  most  common  bridge. 
About  480  out  of  900  arch  bridges  in  Indiana  have  approximately 
this  geometry  (I-b).  Submerged  flow  through  this  type  of  constriction 
was  investigated  for  smooth  and  for  rough  channel  boundaries. 

The  skew  bridges  are  next  in  importance  since  approximately 
35  per  cent  of  all  arch  bridges  in  Indiana  are  skewed.  Among  these 
bridges,  the  most  common  skew  angles  and  the  respective  number  of 
each  structure  in  existance  are:  15°  -  14,  20°  -  80,  and  30°  -  120. 
The  most  common  thickness  factor,  L/b,  is  about  0.50  for  all  bridges. 
Models  of  Geometry  V-b  simulates  the  typical  skewed  arch  bridge 
constriction. 

After  the  one  span  bridges,  two  span  bridges  represent  the  next 
largest  category.  Having  about  100  such  bridges  in  Indiana;  and 
representing  about  11  per  cent  of  all  existing  arch  bridges  in  the 
state.  Models  of  Geometry  VI  simulate  two  equal  and  symmetric  spans 
without  skew  or  rounding. 

Geometry  VII,  represents  two-dimensional  arch  segment  constric- 
tions. Some  model  tests  had  already  been  done  by  a  previous  investi- 
gator, T.  P.  Chang,  but  the  data  taken  had  not  been  analyzed.  These 
experiments  were  complemented  by  additional  testing  and  by  the  comple- 
tion of  the  analysis  of  all  the  data  of  geometry  VH. 


CHAPTER  II 
REVIEW  OF  LITERATURE 

There  is  no  literature  published  to  this  date  which  is  directly 
concerned  with  the  hydraulic  characteristics  of  submerged  arch  bridges. 
Therefore,  papers  on  other  submerged  structures  or  on  the  fundamental 
aspects  of  the  fluid  mechanics  of  contracting  and  expanding  flows 
were  the  only  source  of  literature. 

In  1952  A.  M.  Binnie  presented  a  paper  on  "The  Flow  of  Water 
under  a  Sluice-gate".  By  introducing  the  Froude  number  of  the  flow, 
Binnie  simplifies  relationship  between  conditions  upstream  and  down- 
stream from  a  sluice-gate.  Well-known  expressions  are  derived  for 
alternate  depth  of  flow,  discharge,  and  force  in  the  gate.  A  simple 
explanation  is  presented  for  the  absence  of  waves  on  both  sides  of 
the  gate.  Coefficient  of  contraction  Cc  of  the  issuing  flow  is  given 
as  a  review  of  experimental  and  theoretical  investigations.  Actually, 
Binnie  does  not  mention  Cc,  but  presents  instead  the  gate  opening 
required  to  produce  a  given  Froude  number.  Although  the  paper  pre- 
senting the  problem  uses  the  simplified  one  dimensional  sluice  opening, 

it  is  a  good,  clear  compilation  of  the  present  knowledge  of  the  subject. 

7 
P.  Franke  also  studies  the  one-dimensional  flow  from  a  sluice 

gate  in  his  paper  of  1955  on  "Jet  Contraction  in  Flow  under  Sluice- 
gates in  Rectangular  Channels".  His  analysis  is  based  on  energy 


considerations.  He  computes  the  downstream  depth  in  the  cases  with  and 
without  energy  losses  taken  into  consideration.  Values  of  the  coeffici- 
ent of  contraction  are  given  for  different  ratios  of  upstream  depth  to 
gate  opening.  All  results  are  'expressed  in  non-dimensional  form. 

In  1950,  M.  L.  Albertson,  R.  A.  Jensen,  Y.  B.  Dai  and  H.  Rouse 
published  a  paper  on  "Diffusion  of  Submerged  Jets".  This  paper  deals 
with  the  turbulence  generated  at  the  edge  of  a  free  air  jet  issuing 
from  orifices  and  slots.  Results  are  given  of  measurements  of  the 
distribution  of  the  longitudinal  velocity  and  turbulence  characteristics. 
Of  particular  interest  to  bridge  hydraulics  is  the  discussion  of 

Q 

H.  R.  Henry  (1950)  on  this  paper.  Henry  investigates  the  effect  of 
boundary  conditions  which  are  different  from  those  used  in  the  main 
paper.  These,  for  the  case  of  the  flow  under  a  sluice  gate,  are  the 
effect  of  the  free  surface,  i.e.  gravitational  effect,  and  the  presence 
of  a  solid  boundary,  i.e.  the  flume  floor,  instead  of  a  plane  of 
symmetry  of  a  two-dimensional  jet.  As  a  hydraulic  design  conclusion, 
experimental  discharge  coefficients,  Cd,  for  the  flow  under  a  sluice 
gate  are  given  as  a  function  of  a  dimensionless  gate  opening,  h/B, 
using  a  dimensionless  tailwater  depth,  t/B,  or  the  orifice  Froude 
number  as  a  parameter.  The  experimental  trend  of  the  curves  is  verified 
by  an  approximate  theoretical  analysis. 

In  1957  J*  C.  Gibbings  and  J.  R.  Dixon   presented  a  paper  on 
" Two-Dime nsional  Contracting  Duct  Flow".  Although  this  paper  is 
primarily  for  wind  tunnel  design  dealing  with  flow  through  two-dimensional 
contracting  channels  of  definite  length,  it  presents  a  valuable  potential 
flow  study,  and  a  detailed  method  of  eliminating  the  adverse  velocity 


gradient  along  the  channel  walls. 

Most  of  the  literature  reviewed  so  far  has  dealt  with  the  one- 
dimensional  opening  in  so  far  as  they  use  a  uniform  "infinitely"  wide  • 
slot  as  sluice  gate  opening.  But  the  submerged  arch  bridge  constric- 
tion is  a  two-dimensional  opening  case  where  the  width  of  the  opening 
varies  with  the  height.  This  eliminates  some  of  the  simplifications 
made  in  the  papers  discussed  above,  and  the  present  investigation  has 
to  rely  more  upon  the  experimental  study.  This  is  because  it  is  not 
possible  to  develop  a  simple  mathematical  model  of  the  flow  patterns 
through  two  dimensional  orifices  subjected  to  the  influence  of  gravity. 
However, the  flow  could  in  some  respects  be  analyzed  with  the  help  of 
the  orifice  theories. 

H.  Rouse   discusses  the  characteristics  of  orifice  flow  and  the 

variation  of  orifice  coefficients  with  boundary  geometry. 

12  12 

H.  Rouse   and  A.  H.  Abul-Fetouh   presented  the  paper  "Character- 
istics of  Irrotational  Flow  Through  Axially  Symmetric  Orifices"  in 
1950.  The  authors  conclude  that  use  of  the  relaxation  method  has 
permitted  a  numerical  determination  of  flow  characteristics  to  be  made 
with  sufficient  precision  for  the  problem  to  be  considered  solved, 
though  an  exact  solution  of  the  orifice  problem  has  not  yet  proved 

feasible.  The  coefficient  of  contraction  was  found  to  be  practically 

13 
identical  with  that  evaluated  by  von  Mises  J   for  two-dimensional  flow 

from  slots  over  the  entire  range  of  area  ratio,  and  reasonable  agree- 
ment was  shown  to  exist  between  measurement  and  computation.  Coordi- 
nates of  the  jet  profiles  were  presented  in  tabular  and  graphical  form 
and  were  found  to  differ  appreciably  from  those  previously  adopted 


from  the  two-dimensional  case.  A  composite  dimensionless  chart  was 
also  provided  showing  the  distribution  of  pressure  along  the  boundary, 
center  line,  and  across  the  efflux  section  for  the  various  area  ratios. 

M.  C.  Chaturvedi   a  student  of  H.  Rouse,  presented  in  May  1 963 
a  paper  on  "Flow  Characteristics  of  Axisymmetric  Expansions."  He  has 
studied  the  characteristics  and  dynamics  of  flow  for  four  abrupt 
expansions.  Certain  general  observations  were  done.  Chaturvedi 
observed  that  the  large  scale  properties  of.  motion  are  not  affected 
by  viscous  action.  Thus,  beyond  the  effect  of  viscosity  on  the  condi- 
tions of  stability,  it  has  no  influence  on  the  pattern  of  separation 
and  on  the  turbulence  intensities.  With  regard  to  the  dynamics  of  the 
flow,  it  was  noted  that  the  most  significant  feature  of  the  separation 
is  the  turbulence  production.  Because  of  the  complex  nature  of  the 
phenomena, it  was  not  possible  to  express  the  results  and  the  details 
of  the  flow  characteristics  in  simple  mathematical  terms.  A  simplified 
one-dimensional  analysis  of  the  variation  of  head  loss  was  presented. 

For  the  investigation  of  the  three  dimensional  models  under  sub- 
mergence, the  problems  become  similar  to  those  in  the  investigation  of 
short  culverts.  Here,  several  investigations  have  been  published  and 
good  sources  were  available  in  the  literature  searched. 

The  Highway  Research  Board  Research  Report  15-B  on  Culvert 
Hydraulics,  published  in  1953     »   investigates  the  head-discharge 
relationship  in  short  culverts.     The  existence  of  a  transition  zone 
between  part-full  phase  and  full-flow  phase  was  also  observed  by  the 
investigators. 

Bulletin  126  from  Highway  Research  Board  in  195^       presents  a 


study  on  Culvert-Flow  Characteristicso  Ihe  first  'part  of  the  paper 
covers  a  series  of  typical  water-surface  profiles  which  can  occur  in 
a  culvert;  M.  R.  Carstens  and  A,  R.  Holt  were  the  investigators. 
Ihe  second  part  was  presented  by  R.  E.  Schiller,  Jr.  and  treated  the 
effect  of  various  types  of  inlets  on  flow  through  short  model  circular 

pipe  culverts. 

17 
In  1959t  Hunter  Rouse  and  D.  E.  Ketzler   presented  a  paper 

"Hydraulics  of  Box  Culverts".  The  various  graphs  and  experimental 
evidence  were  based  upon  an  extensive  study  on  box-culvert  character- 
istics. But  the  generality  of  the  approach  in  the  presentation  of 
the  non-dimensional  curves  makes  much  of  the  material  applicable  in 
a  qualitative  sense  to  pipe  culverts  and  submerged  arch  bridge  flow 

as  well. 

1  S 
G.  E.  Myers,  J.  J.  Schauer  and  R.  H.  Eustin   presented  in  March 

1963  a  paper  on  "Plane  Turbulent  Wall  Jet  Flow  Development  and  Friction 
Factor."  They  investigated  the  jet  development,  the  velocity  profiles, 
and  the  wall  shearing  stress  in  a  two-dimensional,  incompressible, 
turbulent  wall  jet.  The  maximum  velocity  decay,  jet  thickness,  and 
shearing  stress  were  predicted  analytically  by  momentum  -  integral 
methods.  Experimental  data  concerning  velocity  profiles,  velocity 
decay,  and  jet  thickness  agreed  well  with  those  of  previous  investi- 
gators. The  wall  shearing  stress  was  measured  by  a  hot-film  technique 
and  the  results  help  to  resolve  a  wide  divergence  between  experimental 
values  of  other  investigators.  All  their  tests  were  conducted  with 
air  blowing  through  a  nozzle. 
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CHAPTER  in 

THEORETICAL  ANALYSIS 

Figures  2  and  3  show  the  centerline  profile  and  the  plan  view 
respectively  of  a  typical  flow  through  a  submerged  arch  bridge  opening. 

In  figure  2  the  approaching  water  surface  is  shown  to  have  a  profile 

19 
of  an  M,  curve  as  defined  by  Ven  Te  Chow  ,  Y  or  Y  is  the  normal 

depth  of  the  unobstructed  channel  flow,  Y.  is  the  depth  at  the  section 

of  maximum  backwater,  Y2  is  the  downstream  depth  taken  in  the  deadwater 

flow  protion  of  the  channel  (see  figure  3)«  The  plan  view  of  the 

channel  constriction  flow  in  figure  3  shows  the  streamlines  of  the 

flow  entering  the  opening.  The  velocity  upstream  of  the  constriction 

is  low  and  the  streamlines  follow  the  walls  very  closely.  The  eddies 

sketched  on  the  upstream  side  of  the  opening  moves  in  a  vertical 

plane  in  a  screw  motion  towards  the  opening.  A  particle  in  the  corner 

would  not  remain  in  this  region  but  due  to  these  eddies  it  moves 

towards  the  opening.  The  eddies  downstream  of  the  constriction,  shown 

on  both  sides  of  the  jet,  are  in  the  horizontal  plane  and  are  less 

active.  Particles  could  often  remain  there  for  longer  periods  of 

time. 
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Dimensional  Considerations 

A  dimensional  analysis  was  used  for  the  design  of  the  experiment 
and  for  the  presentation  and  interpretation  of  the  results.  The 
variables  which  govern  such  an  analysis  are  of  three  kinds:  the  fluid 
properties,  the  kinematic  and  dynamic  variables  and  the  dimensions 
defining  the  constriction  geometry. 
The  variables  are: 

(a)  Fluid  Properties 

1 )  v  ,  the  kinematic  viscosity  of  the  fluid 

2)  p,   density  of  the  fluid 

(b)  Kinematic  and  Dynamic  Flow  Variables 

1 )  g,  acceleration  of  gravity 

2)  Y1?  maximum  water  depth  upstream  of  the  constriction 

3)  Y  ,  the  normal  depth  of  flow  in  the  approach  channel 

4)  V  ,  the  velocity  of  flow  at  normal  depth 

5)  n,  Manning's  roughness  coefficient  of  the  approach 

channel 

6)  Ah,  the  maximum  water  surface  drop  across  the 

constriction 

(Ah=Y1  -  Y2) 

(c)  Properties  of  the  Constriction  Geometry  (see  figures  2, 

4,  and  24  for  definition  of  symbols). 

1 )  A  . ,  the  total  normal  depth  flow  area  at  section  1 • 

2)  A  „  the  normal  depth  flow  area  in  the  opening 

of  the  constriction  (see  Fig.  4).  A.  =  A_  for 
submerged  bridges. 
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3)  L/b,  thickness  factor,  where  L  is  the  length  of 

bridge  in  direction  of  flow,  b  is  the  width  of 
opening  at  the  bottom. 

4)  bL,/A  _,  distance  factor  where  L,  is  the  distance 

d  nZ  d 

between  two  parallel  identical  bridges. 

5)  i>y    ,   wing  wall  angle 

6)  02  ,  skew  angle 

7)  ft  ,   segment  factor  defined  by  ft  =  d/r  where  d  is 

the  distance  between  the  channel  bottom  and  the  center 
of  the  circular  segment. 

8)  N,  number  of  spans 

9)  e,  eccentricity  defined  by  e  =  1  -  c/a,  where  c  and  a 
are  the  width  on  either  sides  of  the  bridge  openings. 

From  the  above  list  of  variables, 

Y1  =  f1  (Yn,  Vn,  n,  Ah,  V.  p,  g,  An|,  A^,  L/b,  -j-S, 

n2 

<f)v   <f)2,   e,p,   N)  (1) 

Buckingham's  theorem  states  that  in  a  physical  problem  including 
n  quantities  in  which  there  are  m  dimensions,  the  quantities  may  be 
arranged  into  (n  -  m)  dimensionless  parameters.  With  the  mass, 
length  and  time  system  of  units,  the  n  -  m  or  14  dimensionless 
parameters  are  as  follows: 


-L.  -  f  UlL      JL2.  n    il   J^   Ah    L 

Yn  ~  2 V  '    v        '  T^  '   Y  2  '  Y  2  '  ^n  '  b  ' 


n 


bL. 

t-5-  ,  tv   02.  e,ft,   N  )  (2) 

*n2       c        ^ 
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V  2  V  Y 

where  s*-  is  the  square  of  the  Froude  number  Fn  and  -&-&•  is  the 
Yng  V 

Reynolds  number  R.  It  is  known  thats 

1 )  Viscous  forces  play  a  negligible  role  in  open  channel 
flow.  So  the  term  R  can  be  neglected. 

2)  The  backwater  superelevation  upstream  is  not  materially 
affected  by  the  shape  of  the  downstream  water  surface, 

A  Vi 

and  the  term  =r—  can  be  neglected. 

n  2 

Furthermore,  by  combining  the  ratios  kjx        and 

A^./Yyj  into  An_/An>  =  M',  and  letting  n  be  constant, 

which  means  that  the  boundary  roughness  in  the  flume  is 

kept  the  same,  then  the  above  equation  simplifies  to: 

Y1   _  .  ^  .,.  L  bLd 


Yn 


f  (Fn,  M\  g.  ~   ,  0V   02,   e,fi,   N  )  (3) 


Of  the  nine  variables  only  two,  the  Froude  number  and  the  contrac- 
tion ratio  describe  the  flow  field.  The  other  seven  variables  in 
the  dimensionless  analysis  describe  the  different  model  geometries 
that  could  be  tested  separately.  For  a  definition  of  the  separate 
geometries  see  section  on  Definition  of  Test  Geometries  in  Chapter  V. 

Definition  of  Channel  Opening  Ratio  and  Channel  Width  Ratio 

In  an  earlier  investigation  on  the  hydraulics  of  arch  bridges 

P.  F.  Biery  defined  the  channel  openings  ratio  (M*).  Referring  to 

figure  4t  the  channel  opening  ratio  is  defined  as  that  portion  of 

the  total  normal  depth  flow  which  can  pass  through  the  bridge  waterway 

without  contraction.  This  is  equivalent  to  the  ratio  A  ./A  <  obtained 

from  the  dimensional  analysis. 

K<  =  q/Q  =  kj^ 
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For  the  submerged  case  where  orifice  flow  occurs,  the  normal 
depth  for  the  flow  is  close  to  or  exceeds  the  height  of  the  model 
openings.  So  for  all  practical  purposes  the  area,  A  ,  of  the  bridge 
opening  can  be  used  instead  of  A  ?t  under  said  conditions.  And  M* 
may  be  defined  for  submerged  bridges  as 
M»  =  q/Q  =  AQ/An1 

Beside  the  Froude  number  (F  )  and  the  channel  opening  ratio  (M'), 
the  channel  width  ratio  (M)  is  perhaps  the  most  significant  variable 
in  the  problem.  The  channel  width  ratio  is  defined  as  the  ratio  of 
the  model  opening  width  at  the  bottom,  b,  to  the  width  of  the  channel,  B» 
M  =  b/B 

Both  M  and  M1  are  ratios  related  to  the  geometry  of  the  constric- 
tions, while  the  Froude  number  represents  the  effect  of  gravity  upon 
the  state  of  flow. 

Orifice  Flow 

A  distinction  between  the  occurrence  of  free  surface  flow  and 
orifice  flow  is  essential,  and  some  means  is  needed  for  predicting 
when  the  different  types  of  flows  occur. 

Orifice  flow  will  occur  when  the  head  over  the  crown  of  the  arch 
exceeds  the  velocity  head  plus  the  head  loss  in  the  orifice  type 
opening. 

Putting  this  statement  into  an  equation  gives,  for  the  limiting 

case:     ,  0 

V  V 

■ir  +  ( 72  - 1  >  ir  ■  ^  - r  w 
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where  (  — s-  -  1  )  •** —   is  the  orifice  head  loss. 

C  g 

v 

D»n  —   ir  =  Yi  - r 

V 

and  from  the  continuity  equation: 

V1   BY1  =  V3  A3  =  V3  °c  Ao 

where  A  is  the  opening  area  of  the  model. 

0  2       2       2 

1        V     B     Y1 
Then       -~    —         ?     ,     =     Yi   -  r  (5) 

C         2g       C    A 

V  &  CO 

V1  r  2       A       2 

rearranging   (5):     2p  v      =     (  1   -  —  )     C^  (  —  ) 

Hence         F.2  =  2  C,2   (   1    -  ~  )      (  Ao   )2  (6) 

and  there  will  be  submergence  when: 

"i2<  2  cd2  < 1  -  ^  >  <  if  >2  <7) 

An  average  experimental  value  of  C,  is  0.57.  (See  Chapter  VII) 

The  Equation  of  Discharge  for  Orifice  Flow 

With  reference  to  Figure  5  and  neglecting  the  velocity  of  approach, 
the  discharge  was  calculated  as  follows: 


w 


VdA  (8) 


where 


V  =  Cd  V  2g  (Y1  -  h) 
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h  =  r 

and  A  =  2   J     V    r2  -  h2        dh 

'     h  n  6 
r 

hence       Q     =  2      /     V    Cd2     2g   (  Y1   -  h)    V  r2  -  h2         dh 
0^  7 


dh 


=  2  Cd  V2g         J     (T1    -  h)^  (   r2  -  h2   )* 


An  approximate  solution  to  the  above  integral  is  found  by  expand- 
ind  each  term  of  the  integrand  into  a  binomial  series,  multiplying 
the  two  series  and  integrating  term  by  term.      The  result,   using  the 
5  leading  terras,   is 


Q     =  0.^019  C,V2g^     Y*  b2  [1    -  0.2136   (  7-  )   -  0.0321 6   (  £-  )2 
a  i1  i1 

-0.0112       (  $-  )3     -     0.0053^  (  f-  f 1  (10) 

l1  i1  J 


or  Q     =  G,   Y^2     b     T 


where 


C1    =  0.401 9 V2g       Cd  =  3.22  Cd 


and 


T  =  1    -  0.2136    (  $-  )  -  0.03216   (  f-  )2  -  0.0112   (  $-  ) 
I1  .  i1  i1 

-  0.0053^   (  7-  )4  ••■ 

The  above  equation  is  valid  as  long  as  the  constriction  is  sub- 
merged, which  means  for  Y,>r. 
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Backwater  Ratio  Equation  for  Submerged  Inlet  with  Orifice  Flow 

The  backwater  ratio  is  defined  as  the  ratio  of  the  maximum  center- 
line  depth  to  the  normal  depth  of  flow.  See  Figure  2. 

An  expression  for  the  backwater  ratio  may  be  obtained  from  the 
theoretical  discharge  equation  for  orifice  flow.  The  equation  for 
the  discharge  through  the  constriction  was  found  to  be: 

Q  =  G,  Y^  b2  T  (11) 

The  discharge  in  the  approach  channel  is: 

Q  =  V  A  =  F  A/g~  BY3'2  (12) 

^    n  n    n  v  6     n  v   ' 

Equating  the  two  discharge  equations  11  and  12  : 

C1  It*  b2  T  =  Fn<v/g"   Yn3/2  B  (13) 


Sind  Yn     1    K     BYn         Y„   1      R       Yn  1 


but 


b  *  b    k  A     k  M» 

o 


where   k  =  — 


Now  rearranging  the  equation  (14)  : 

i  -  <  rw  F„  •  -h;  ■  +  >2  <•« 

n 
or  in  general  a  relation  of  the  following  form  would  be  expected 
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f  =  c  (  1  )2 

n 


(16) 


It  has  been  found  empirically  that,  instead  of  equation  (16),  the 
exoerlmental  results  can  conveniently  be  represented  by  a  function  of 
the  form: 

=  y~   ~  1    -    f  (  vT  ) 

Y      n 

n 

where  h.   is  the  backwater  superelevation. 

Experimental  results  for  the  orifice  flow  have  yielded  the  following 
relation  (see  Chapter  VII) 


Y« 


=  1  +  1.18 


tkf 


0.90 


(1?) 


J-  =  1  +  1.18  (  r£  ) 


n  >1.80 


(18) 


The  experimental  result  for  the  backwater  ratio  equation  for  free  surface 
was  developed  by  P.  F.  ;;iery  and  found  to  yield  the  relation 

Y 


Y~   =  1  +  c.^7 


(     —  V 


J.  39 


(19) 


-i.  =   1  +  0.^7  (  7TT  ) 


n  x2.26 


(20) 


Dimensionless  Equation  for  Orifice  Type  Flow 


When  the  approach  velocity  to  the  constriction  is  neglected, 
orifice  type  flow  is  a  good  approximation  to  the  flow  condition 
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existing  through  a  fully  submerged  arch  bridge.  In  all  the  tests 
performed  on  the  submerged  bridge  openings  the  approach  velocity  was 
very  low  relative  to  the  velocity  in  the  discharge  jet,  and  the 
assumption  appears  to  be  safe.  The  rate  of  flow  through  a  submerged 
semi-circle  type  bridge  opening  can  thus  be  written  as  follows  (ref. 
Figures  2  and  5): 


where 


and 


•Q  =  discharge  through  the  opening 
d  =  diameter  of  the  opening 

L.  =  the  depth  of  the  backwater 


C ,  =  coefficient  of  discharge  for  type  of  inlet 

Mow  rearranging  equation  (21): 

2   2 
p     r 

g  d 

By  dividing  by  d,  and  taking  the  square  root  of  both  sides  of  (22), 
the  following  dimensionless  equation  is  obtained 


t9™  ■  ^-    <  r  >*  ™ 


Cd7T 
where  — — —      is  a  constant. 

732" 
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The  previous  equation  may  be  written  in  the  general  form 


Y1 

372         r  i  T 


g3  d- 


Vlhich  is  used  in  plotting  experimental  results 

Head  Loss  Due  to  Bridge  Constriction 

An  equation  for  the  head  loss  due  to  the  constriction  effect  may- 
be derived  from  an  energy  equation  written  between  the  section,  1 , 
of  maximum  depth  and  section  4  where  the  flow  has  regained  its  normal 
depth.  With  reference  to  Figure  2. 

V  V 

S0L1J*  +  T1+Q!1     "2i    =  Y4+Q4     ~    +  E1-4  <28> 

where       S     =  slope  of  channel  bottom 
o 

L-ik  =  distance  between  sections  1    and  k 

Y..    and  Y^  =  water  depth  of  sections  1   and  4  respectively 

V.   and  VY  =  mean  velocity  at  sections  1    and  4 

GL    andQj.  are  the  respective  kinetic  energy  coefficients 

hence       S     L.    ,,  =  the  difference  in  bottom  elevation  between  the  two 
sections 

and  E1    k  =  the  total  head  loss  between  sections  1    and  4. 

The  energy  loss  between  sections  1    and  4  may  be  considered  to  be 
made  up  of  two  parts: 

(a)     The  head  loss  that  would  exist  if  uniform  flow  existed 
between  these  sections,  which  would  be  equal  to  the 
change  in  elevation  between  these  two  sections  as  in 


2k 


uniform  flow  the  gravity  forces  exactly  balance  the 
resistance  forces, 
(b)  The  additional  head  loss  due  to  the  non-uniformity  of 

the  flows  namely  the  combined  effect  of  the  contraction, 
separation  and  expansion  of  the  flow  as  it  passes 
through  the  constriction.  The  total  head  loss  may  be 
written  as: 

E1-4  =  hf1    +hf2  (29) 

where  hf1    is  the  head  loss  in  the  unobstructed  channel  and  hence 
equation  29  becomes 

E1 J+  =  So  h-Ur  +  hf2 

It  has  been  suggested  by  C.  Izzard  £k)that  h»  can  be  expressed  as  a 

function  of  the  velocity  head  in  the  opening 

V2 
hf=K-fj  (31) 

where 

V  =  -S- 
o    A 
o 

K  is  the  additional  head  loss  coefficient  due  to  the  bridge 

constriction.  Now  equation  28  simplifies  to: 

2  2        2 

V.         V.       V 

Y1  -Y^=  G4  Tz     "  al  4g  +K  "fe  (32) 


setting     Y1   -  1u  -  h1 *     and  rewritten  equation  32  gives 

V   2  r  Vk      9  V1       9  i  V   2 

V  =  «-!i  +[a*<?t)2-ai  Cf  ) 2]      -fc  03) 


QSjl  and  Q!.  can  be  assumed  to  be  equal,  and  by  using  the  relationship 
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between  the  velocities  and  areas, equation.  33  can  be  written: 


V 


(  sf  )2  -  (  f  )2 

A4      A1 


V 


o 

2g 


(3*0 


Divide  by  -s—  and  solve  equation  3^  w.r.t.  K  gives  an  equation  for 
*-g  ' 

the  additional  head  loss  coefficient  due  to  the  bridge  constriction, 
for  submerged  conditions. 


K  = 


V 
V02/2g 


-a 


(-^  )2  -  (  f  )2 

A4       A1 


(35) 


A  A 

The  ratio  t—  is  the  same  as  r—  (ref.  figure  2)  and  hence  equal 

U-  n 

to  the  channel  opening  ratio  M*. 
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CHAPTER  IV 
EXPERIMENTAL  APPARATUS 

Flune 

All  tests  were  performed  in  a  tilting  flume  5  feet  wide,  2  feet 
deep  and  64  feet  long.  The  flume  was  made  of  welded  steel  plates 
coated  with  an  epoxy  resin  paint.  Six  screw  jacks,  three  on  each  side, 
placed  underneath  the  flume  made  it  possible  to  set  the  flume  bed  to 
any  desired  slope  with  an  accuracy  of  t   0.000025  feet/feet.  £ach 
pair  of  jacks  had  a  different  gear  ratio  and  was  placed  in  such  a  way 
that  96  turns  of  a  common  drive  shaft,  connected  to  an  electric  motor, 
would  raise  the  flume  1,  2,  and  3  inches  respectively,  the  upstream 
pair  of  jacks  raising  1  inch.  The  electric  motor  was  operated  by  three 
switch  buttons  -  raise,  lower  and  stop.  A  revolution  counter,  calibrated 
to  indicate  the  slope,  made  it  possible  to  adjust  the  flume  to  the 
proper  slope  within  minutes.  Figure  6  shows  a  plan  view  of  the  jacks 
and  gears  and  Figure  7  shows  jack  detail. 

At  the  upstream  end  the  flume  was  hinged  to  a  forebay  consisting 
of  an  S  foot  by  10  foot  headbox  lined  with  sheet  metal.   To  assure  a 
smooth  flow  condition  from  the  headbox  into  the  flume,  the  box  was 
equipped  with  a  stone  baffle,  several  screens  and  a  smooth  elliptical 
transition  into  the  flume. 

At  the  discharge  end  of  the  flume  an  adjustable  weir  type  tail- 
gate provided  control  of  the  depth  of  the  flow.  It  could  be  raised 
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or  lowered  by  means  of  a  calibrated  wheel  (k   turns  per  inch). 
Figure  8  shows  tailgate  constructions  A  catchment  box  eliminated  any 
splash  as  the  flume  discharged  the  water  over  the  tailgate  and  into 
the  sump.  For  a  more  complete  description  of  the  testing  flume  see 
reference  2.  Figure  9  shows  apparatus  arrangement. 

Water  Supply  and  Measurement  System 

A  recirculating  sump  provided  the  water  to  the  flume  by  means  of 
two  centrifugal  pumps  of  capacity  2000  GPM  and  300  GFM.  These  fed 
the  headbox  through  a  6  inch  and  a  3  inch  pipeline.  The  flow  was 
measured  by  two  venturi  meters  connected  to  a  30  inch  manometer  filled 
with  mercury,  (Sp.  Gr.  13.6)  for  the  6  inch  line  from  the  2000  GPM 
pump,  and  to  a  60  inch  manometer  filled  with  tetrabromoethane  (Sp. 
Gr.  2.95),  for  the  3  inch  line  and  300  GPM  pump. 

The  calibration  of  the  venturi  meters  were  done  by  means  of  a 
concrete  side  channel  and  a  large  weighing  tank.  The  flow  through 
the  pipelines  was  controlled  by  gate  valves.  Calibration  charts 
were  made  for  both  venturi  meters.  (Figures  10  and  11  calibration 
charts  for  6"  and  3"  pipes).  The  manometers  were  placed  close  to  the 
flow  controlling  valves  so  their  deflections  were  easily  visible  when 
adjusting  the  valves. 

To  eliminate  undesirable  surface  waves  which  reflected  back  and 
forth  in  the  upper  part  of  the  flume,  a  skimming  board  was  placed  at 
the  entrance  in  the  flume. 
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FIG    7   JACK      DETAIL 


FIG    8    TAIL     GATE     CONSTRUCTION 
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Test  Station 

To  be  able  to  change  models  while  the  water  was  running,  to  put 
the  models  into  the  flow  at  the  particular  station  -  at  a  particular 
angle  -  was  essential  in  order  to  secure  uniform  conditions  at  the 
testing  section  in  the  flume.  This  problem  was  solved  by  the  use  of 
two  1 M  x  1 "  x  1 8"  angle  irons.  They  were  bolted  to  the  sidewalls  of 
the ' flume  in  such  a  way  that  they  easily  could  be  removed  when  calibra- 
tion of  the  normal  flow  depth  had  to  be  done.  Figure  12  shows  the  test 
section  with  a  model  in  place.  Sach  angle  iron  had  two  holes  where 
bolts  leading  through  the  model  could  be  fastened. 

To  change  models,  the  four  bolts  holding  the  model  were  loosened 
and  the  model  was  lifted  out.  The  new  model  was  slid  down  behind 
the  angle  irons  into  a  perfect  position  at  the  test  station  and  the 
bolts  were  tightened  to  hold  it  in  place. 

Measuring;  Equipment 

Instrument  Carriage 
An  aluminum  instrument  carriage,  mounted  on  two  stainless  steel 
guide  rails  running  the  length  of  the  flume,  was  installed  so  that 
the  bottom  of  the  flume  could  serve  as  a  reference  plane.  The  four 
wheels  of  the  carriage  were  grooved  to  give  a  precision  fit  on  the 
rails.  One  of  the  wheels  had  a  set  screw  to  keep  the  carriage  from 
moving  when  an  experiment  was  conducted.  A  surveyors  tape  was  mounted 
along  the  whole  length  of  the  flume  and  an  indicating  point  on  the 
carriage  served  as  a  reference  point  for  all  longitudinal  measurements. 
Another  surveyors  tape  was  installed  in  the  transverse  direction  of 
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FIG.     12a        TEST       SECTION 


FIG.    12b       TEST     SECTION     DETAIL 
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the  flume  on  a  grooved  rail  of  the  instrument  carriage.  This  rail 
served  as  a  guide  for  the  slide  on  which  a  point  gage  and  a  Prandtl 
tube  were  mounted.  This  slide  also  had  an  indicating  point  serving 
as  reference  point  for  transverse  measurements.  3y  means  of  these 
two  tapes  any  point  in  the  flume  could  be  located  in  the  horizontal 
plane.  The  instrument  carriage  was  also  equipped  with  a  florescent 
lamp,  several  electrical  outlets  and  a  wide  "}/k   inch  plywood  strip  on 
one  side  to  serve  as  a  desk  for  notekeeping  etc.  The  carriage  could 
easily  be  moved  to  any  location  along  the  flume  and  locked  in  a 
specific  position  with  an  accuracy  of  0.005  ft.  (Figure  13  shows  a 
top  view  of  the  carriage.) 

Point  Gage 

All  water  surface  elevation  measurements  were  obtained  by  means 
of  an  electric  indicating  point  gage.  It  was  mounted  on  the  instrument 
carriage  slide.  This  could  be  positioned  with  an  accuracy  of  +  0.002  ft, 
in  the  transverse  direction  of  the  flume.  The  staff  of  the  pointgage 
was  marked  in  milli-neters  and  was  equipped  with  a  vernier  which 
allowed  a  reading  accuracy  of  t   0.05  ™i  in  the  vertical  direction. 
The  accuracy  of  the  water  surface  elevation  measurements  varied 
between  0.1  mm  to  2.0  mm  depending  on  the  smoothness  of  the  free 
surface. 

The  procedure  used  when  the  free  surface  was  not  smooth  was  to 
observe  the  time  intervals  of  the  blinking  light  of  the  point  gage. 
The  light  is  on  when  the  measuring  point  is  at  or  below  the  water 
surface.  When  the  light  was  blinking  on  and  off  with  equal  time 
intervals,  that  is  the  measuring  point  was  in  and  out  of  the  water 
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FIG.    13       TOP     VIEW     OF       INSTRUMENT       CARRIAGE 


FIG      14         POINT      GAGE      AND      PRANDTL     TUBE 
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for  an  equal  duration,  the  reading  of  the  vernier  was  recorded. 

Prandtl  Tube 

The  flow  velocity  measurements  were'  obtained  by  means  of  a  1  /3ff 
O.D.  Prandtl  tube  of  standard  design.  The  Prandtl  tube  was  mounted  on 
a  vertical  staff  which  was  marked  in  0.2  cm  and  had  a  vernier  allowing 
a  vertical  positioning  within  1  0.01  cm  of  accuracy.  A  set  screw 
prevented  the  position  to  change  during  test  measurements.   The 
support  of  the  Prandtl'  tube  permitted,  in  addition  to  vertical  and 
traverse  motions,  a  rotation  so  that  the  horizontal  portion  of  the 
probe  could  be  alligned  according  to  the  direction  of  the  flow  indicated 
by  a  freely  swinging  stiffened  tuft  attached  to  the  tube.   Figure  14 
shows  point  gage  and  Prandtl  tube.  The  amount  of  rotation  was  measured 
as  an  angle  by  means  of  an  indicator  sliding  along  a  fixed  protractor. 

Since  the  tip  of  the  Prandtl  tube  could  be  aligned  with  the 
streamline  direction  in  the  horizontal  plane  only,  the  question  arose 
as  to  the  reliability  of  the  measured  velocity  whenever  it  had  a 
vertical  component.  In  order  to  establish  the  directional  sensiti- 
vity of  the  prandtl  tube,  a  series  of  velocity  measurements  were  taken 
varying  the  angle, y,   between  tip  of  the  Prandtl  tube  and  the  direction 
of  a  uniform  flow.  It  was  established  that  the  velocity  reading  was 
within  about  jt   accuracy  when  the  tip  of  the  Prandtl  tube  was  directed 
within  17°  of  the  flow  direction.  Figure  1^  shows  graph  of  directional 
sensitivity  of  Prandtl  tube. 

Variable-Reluctance  Pressure  Transducers 
The  small  magnitude  of  the  pressure  differentials  to  be  measured 
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with  the  Pradtl  tube  made  it  necessary  to  employ  a  Variable-Reluctance 
Differential  Pressure  Transducer  (model  P7D,  Pace  Eng.  Comp. )  instead 
of  the  conventional  liquid  manometer.  In  the  low  velocity  regions  a 
pressure  transducer  having  a  range  of  t   0.1  psid.  was  used.  In  ihe 
higher  velocity  regions  -  as  in  the  discharge  jet  of  the  models  -  a 
range  at  1-0.5  psid.  was  needed. 

The  range  of  the  pressure  transducer  was  varied  by  changing 
diaphragm  (membrane)  in  the  pressure  cell,  Figure  16  shows  sketch 
of  velocity  pressure  transducer  system.  The  thickness  of  the  diaphragm 
governs  the  pressure  range;  thicker  diaphragm  are  used  for  larger 
differential  pressure  ranges.  The  leads  from  the  static  and  stagnation 
points  on  the  Prandtl  tube,  filled  with  Meriam  No.  D-2930  fluid  of 
controlled  specific  weight  (sp*  gr.  1.00)  and  low  surface  tension, 
were  connected  to  each  side  of  the  transducer  diaphragm.  In  order  to 
prevent  fluctuation  of  the  flow  to  influence  the  pressure  differential 
and  hence  the  voltmeter  reading,  the  leads  contained  coiled  lengths 
of  tubing  -  giving  reservoir  action.  Shorter  pieces  of  leads  with 
much  smaller  diameter  were  also  introduced  to  help  dampen  the  pressure 
fluctuations. 

The  principle  of  operation  of  the  diaphragm  type  variable- 
reluctance  pressure  transducer  is  as  follows: 

A  diaphragm  of  magnetically-permeable  material,  supported  between 
two  symmetrical  E  core  indutance  assemblies,  Figure  1?,  completed  a 
magnetic  circuit  with  each  E  core.  The  diagram  deflected  when  there  was 
a  difference  in  pressure  between  the  two  input  lines.  This  increased 
the  gap  in  the  magnetic  flux  path  of  one  core,  and  decreased  the  gap 
equally  in  the  other. 

The  magnetic  reluctance  varied  with  the  gap,  determining  the 
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inductance  value,  so  that  the  over-all  effect  of  diaphragm  motion 

resulted  in  a  change  in  inductance  of  two  coils  L1  and  L?.  The 

inductance  ratio  L,  /   was  conveniently  measured  in  a  bridge  circuit 

L2 
which  produced  an  output  voltage  proportional  to  the  pressure  difference, 

Between  the  Pressure  Transducer  and  the  voltmeter  a  Carrier- 
Demodulator  (model  CD10,  Pace  Eng.  Comp. )  provided  a  carrier  excita- 
tion of  3000  cps.  (Figure  18  Carrier-Demodulator).  It  was  operated 
from  the  115V60Cy  line.  The  span  control  knob  provided  a  scale  adjust- 
ment and  the  zero  control  knob  set  the  voltage  output  corresponding 
to  zero  pressure  difference. 

A  Junior  Volt  Ohmyst  Voltmeter  (model  WV-77A,  RCA)  was  used. 
This  also  had  zero  setting  and  five  scale  adjustments  besides  a  +  and 
-  volt  selector  switch.  The  voltmeter  was  placed  on  a  rotating  disk 
to  obtain  a  flexible  position  for  the  reader,  who  had  to  read  the 
instrument  from  several  different  positions  during  calibrations  and 
testing. 

Three  diaphragms  were  available  which  could  be  used  for  pressure 
differential  measurements  up  to  0.1,  0.5  and  1.0  psi.  A  diagram  was 
prepared  for  the  selection  of  the  proper  diaphragm  for  a  given 
velocity  range.  The  velocity  is  related  to  the  pressure  differential 
by  the  equation 


V  =V2gAp/r 
where  A  p  is  the  pressure  differential  measured  between  the  dynamic 
and  static  openings  of  the  Prandtl  tube.  The  pressure  differential 
was  plotted  vs.  the  velocity  on  Figure  ^9,   where  the  range  of  the 
several  diaphragms  were  indicated. 
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-OUTPUT   CONNECTOR: 


TRANSDUCER    INPUT  CONNECTOR: 


<* 


POWER    INPUT: 
115  V      60~ 


FIG     18     CARRIER    -    DEMODULATOR 
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Calibration  of  Velocity  Probe 
By  using  the  relation  between  the  dynamic  and  static  head 

r—  =  h  (where  V  is  velocity  in  fps.,  h  is  feet  of  water  and  g  is  the 
acceleration  of  gravity,)  a  static  calibration  could  be  performed. 

Two  reservoirs  were  connected  to  the  stagnation  and'  static 
openings  of  the  Prandtl  tube  respectively.  The  reservoirs  were 
connected  through  a  valve.  A  datum  level  was  measured  by  an  electric 
point  gage  when  the  connecting  valve  was  open,  and  the  voltmeter 
needle  was  adjusted  to  the  zero  mark.  The  pressure  was  then  equal 
on  both  sides  of  the  diaphragm  in  the  pressure  transducer.  The  valve 
connecting  the  reservoirs  was  closed  and  distilled  water  was  added 
to  increase  the  level  in  the  reservoir  connected  to  the  stagnation 
opening.  The  water  level  was  recorded  for  even  increments  in  the 
voltmeter  reading.  Figures 20  and  21  show  the  calibration  setup. 
By  means  of  a  least  square  subroutine  of  an  electronic  computer 
(see  Appendix  C)  the  relation  between  the  output  voltage  and  velocity 
was  determined  from  the  measured  pairs.  This  realtion  was  found  to 
be  a  constant,  C,  and  the  velocity  was  found  from  the  expression 
V  (fps.)  =  C^/volt. 

A  new  calibration  test  had  to  be  made  each  time  the  diaphragm 
was  changed.  Repeated  calibrations  without  changing  the  diaphragm 
gave  the  same  calibration  curve.  Typical  calibration  curves  for 
probe  are  shown  in  Figure  22  where  velocity  versus  voltage  are  plotted. 
During  the  experiments  it  was  necessary  to  measure  both  high  and  low 
velocities.  To  omit  the  timeconsuming  process  of  calibration,  two 
individual  pressure  transducers  were  used  with  a  range  of  0.5  psid. 
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FIG     20        CALIBRATION      APPARATUS      FOR     PRANDTL      TUBE 


47 


FIG.     21        DETAIL      OF      CALIBRATION  APPARATUS 
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Fig  22    Typical     Calibration     Curves    for    Probe 
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and  0.1  psid.  respectively.  They  were  mounted  on  a  lucite  plate  in 
such  a  way  that  one  could  switch  from  one  to  the  other  in  a  matter 
of  seconds.  The  lucite  plate  was  mounted  on  the  slide  of  the 
instrument  carriage  as  shown  in  Figure  23,  which  also  shows  the  change 
in  position  of  the  valves  as  one  or  the  other  pressure  transducer 
was  in  use. 

Procedure  for  Measuring  of  the  Velocity 
After  checking  that  no  air  was  trapped  in  the  pressure  leads, 
equal  pressure  was  applied  on  both  sides  of  the  diaphragm  by  connecting 
both  leads  to  the  same  reservoir.  Then  the  voltmeter  needle  was 
adjusted  to  zero  on  the  proper  scale  needed  for  the  expected  velocity 
range  in  the  experiment.  The  valves  to  the  equal  pressure  reservoir 
were  closed,  the  Prandtl  tube  was  inserted  into  the  flow  and  the  cap 
covering  the  stagnation  and  static  openings  was  removed.  Now  the 
valves  connecting  the  stagnation  and  static  pressure  tubes  with  the 
pressure  transducer  were  opened  to  introduce  the  pressure  difference 
upon  the  diaphragm  and  thus  make  a  change  in  the  output  voltage. 
The  Prandtl  tube  was  placed  in  the  desired  position  and  the  voltmeter 
reading  was  recorded. 

After  the  test  measurements  were  taken,  the  instruments  were 
again  checked  against  the  equal  pressure  reservoir  to  see  that  the 
voltmeter  read  zero  for  the  same  pressure  on  both  sides  of  the 
diaphragm. 

The  calibration  apparatus  and  the  technique  of  the  use  of  the 
pressure  transducer  for  velocity  measurement  were  developed  by 
A.  A.  Sooky  as  part  of  another  research  project. 
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RIGHT     CELL     CONNECTED     FOR      HIGH       VELOCITY      MEASUREMENT 


LEFT      CELL      CONNECTED      FOR      LOW      VELOCITY      MEASUREMENT 
FIG.    23  PRESSURE         CELLS 
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CHAPTER  V 

DEFINITION  OF  TEST  GEOMETRIES  -  SELECTION  OF  TESTS 

Nine  independent  variables  were  considered  in  the  dimensional 
analysis  to  determine  the  dependent  variable  Y.  /Y  .  Of  these  nine 
variables,  seven  describe  different  types  of  constriction  geometries, 
the  other  two,  the  Froude  number  and  the  contraction  ratio  describe 
the  flow  field  and  the  amount  of  contraction.  These  seven  geometric 
variables  were  used  to  define  seven  types  of  geometries  that  were 
tested  separately.  These  types  of  geometries  are  defined  below. 


Geometry  I-a 
Two-Dimensional  Semicircle  Arch  Bridge  Constrictions: 
(Figure  24)  The  characteristics  of  this  type  of  geometry  are 
the  following: 

a)  Since  the  model  is  two-dimensional,  r  =  0 

o 

bLd 

Single  bridge  case,  L ,  =  0  and  t —  =  0 
a        n2 

No  wingwalls,  according  to  definition  of  <f>. ,  <f>.,   =  90° 
None-skew  case,  according  to  definition  of  0_,  ?L,  =  0 
Semicircular  case,  d  =  0  and  j3  =    7=0 
One  span  case,  N  =  1 
g)  No  eccentricity,  e  =  0 
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V 
h)  Froude  number,  F„  =     1  designated  as  an  independent 

variable. 


i)  Constriction  opening  ratio  M1,  designated  as  an  independent 
channel  opening  ratio  M* 

L   bLd 
Hence,  seven  independent  descriptive  variables,  —  ,  r—  ,  e  , 

D         n2 

Y1 

$.,   0o./5  and  K  in  ^e  dimensionless  equation:      (3)..  s-    =  f  (F  , 

n 

L   DLd 
M*,  r  ,  7 —  ,  $*,  $2*   e»/5*  ^   are  designated  as  constants  and  we 
n2 

get  the  relation: 

YL=  f  (F  .  M») 
n 

Geometry  I-b 

Three-Dimensional  Semicircular  Arch  Bridges: 


This  geometry  differs  from  the  previous  one  only  in  that  L  is 
varied,  thus  the  parameter  r  describing  the  length  of  the  model  in 
the  direction  of  the  flow  is  a  variable  and  the  dimensionless  equation 
(3)  for  this  case  becomes 

£-f  CF„.  M'.  i) 

n 

Geometry  II 
Dual  Parallel  Three-Dimensional  Arch  Bridge  Constrictions: 
This  geometry  consists  of  two  identical  bridges  of  geometry  I-b, 
placed  at  a  distance  L,  apart,  measured  center  to  center.  One  new 
variable,  L,,  is  introduced,  which  is  characterized  by  the  parameter 
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BL, 

t—  (Note:  For  submerged  bridge  constrictions,  A  2  =  A,.,  see  Chapter 
n2 

III).  Equation  (3)  simplifies  for  this  case  to: 

Y1  T   bLH 

n  n2 

However,  in  order  to  simplify  the  testing,  the  length  ratio  — 

was  limited  to  the  case  of  r-  ='  0.25i  so  that  this  quantity  was  a 

constant  for  the  tests,  and  equation  (3)  reduces  to: 

Y,  bL 

^=f  (Fn.  M«.jJ) 

n  nZ 


Geometry  III 
Three-Dimensional  Arch  Bridge  Constriction  with  Wingwalls: 

The  geometric  characteristics  were  as  follows: 

a)  £=  0.25  e)  /5=  0 
bL 

b)  r-a  =  o  f )  n  =  1 
V 

c)  01  is  a  variable  ^  =  90°,  60°,  hf,   or  30°      g)  e  =  0 

d)  0-  =  0°  h)  F  is  a  vari- 


able 


i)  M*  is  a  variable 


Hence,  equation  (3)  can  be  simplified  to: 

Y1 

I 

n 


Y1 

Yi=  f  (Fn.  M\  0t) 
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Geometry  IV 
Two-Dimensional  Semicircular  Arch  Bridge  Constrictions  with  Eccentricity; 
The  geometric  characteristics  were  as  follows: 

a)  —  =  0  d)  (L  =  0  g)  e  is  a  variable 

bL 

b)  '--£  =0  e)  B  =  0  h)  F  is  a  variable 

V  n 

c)  01  =  90°  f)  N  =  1  i)  M1  is  a  variable 

And  the  equation  (3)  can  be  reduced  to: 
Y 

— 1  =  f  (y       w'   e  ) 

n 

Geometry  V-a 
Two-Dimensional  Semicircular  Arch  Bridge  Constrictions  ~..:ith  Skew: 
The  geometric  characteristics  were  as  follows: 

a)  r  =  0  d)  $~   is  a  variable     g)  e  =  0 

bL. 

b)  T-2-  =0  e)  j3  =  0  h)  F  is  a  variable 

n2 

c)  L   =  90°  f )  N  =  1  i)  M'  is  a  variable 

And  the  equation  (3)  can  be  reduced  to: 

Y, 

YL=  f  (Fn,  !'•,  <f>2) 
n 

Geometry  V-b 

Three -Dimensional  Semicircular  Arch  Bridge  Constrictions  With  Skew: 

This  geometry  is  the  same  as  the  previous  one,  except  that  the 

length  of  the  constriction  is  allowed  to  vary,  the  parameter  -  is  thus 

a  variable  and  the  dimensionless  equation  (3)  can  be  simplified  to: 

Y    1  un'  r<  »  v2*   b  ' 

n 
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Geometry  VI 
Three-Dimensional  Two-Span  Semicircular  Arch  Bridge  Constrictions: 
The  geometric  characteristics  were  as  follows: 

a)  r  =  0.50  f)  N  =  2,  Two-span  case 

b)  r-S  =  0  g)  e  =  0 

n2 

c)  0^  '=  90°  h)  Fn  is  a  variable 

d)  02  =  0°  ,      i)  M'  is  a  variable 

e)  P   =  0  °;  ?"  10 

The  equation  (3)  can  be  reduced  to: 

Y    r  v*nt  n  ; 

n 

Geometry  VTI 
Tuo-Dinensional  Segment  Arch  Bridge  Constrictions: 
The  geometric  characteristics  were  as  follows: 

a)  r  =  0  e)  Q  is  a  variable 

bL, 

b)  j-2-  =  0  f )  N  =  1 

n2 

c)  ^  =90°  g)  e  =  0 

d)  02  =  0°  h)  Fn  is  a  variable 

i)  M*  is  a  variable 
Hence,  equation  (3)  reduces  to: 

n 
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Submerged  Test  Geometries: 

The  test  geometries  defined  and  sketched  in  the  previous  section 
give  an  overall  picture  of  the  variety  of  models  that  were  considered 
in  the  arch  bridge  testing  program.  Only  five  of  these  geometries 
were  tested  for  submergence.  (See  Introduction,  for  background  of 
model  selection) 

Geometry  I-a  was  selected  to  be  tested  first.  This  two-dimensional 
geometry  was  the  simplest  possible  arch  model.  It  was  tested  in  the 
smooth  channel  to  give  some  basic  information  on  the  velocity  dis- 
tribution and  on  the  coefficients  of  contraction,  of  velocity  and 
of  discharge.  It  also  served  as  a  guide  for  further  investigation 
on  the  more  complicated  constrictions. 

The  three  dimensional  geometry  I-b  was  the  next  step.  These 
models  were  also  tested  in  the  smooth  boundary  channel  and  the 
results  are  presented  in  non-dimensional  form.  After  the  exploration 
of  the  test  range  was  completed,  and  a  suitable  way  of  analyzing  the 
data  was  found,  the  flume  boundaries  were  changed  from  smooth  to 
rough. 

Roughness 

The  rough  boundary  was  obtained  by  placing  £  inch  aluminum  rods 
along  the  flume  bottom  and  side  walls.  Along  the  bottom  a  layer  of 
longitudinal  bars  were  placed  12  inches  on  center  and  a  top  layer  of 
transverse  bars  6  inches  on  center.  Along  the  side  walls  one  layer 
of  vertical  bars  6  inches  on  center  was  placed  £  inch  from  the  wall . 
The  bottom  layer  of  bars  were  tied  together  with  wire.  The  vertical 
bars  were  tied  at  the  bottom  to  the  transverse  bars  and  clamped  to 
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the  walls  above  the  free  surface.  Figure  1  shows  the  roughness  bars 
placed  in  the  flume.  The  roughness  was  calibrated  and  Mannings 
roughness  coefficient,  n,  was  found.  The  average  Mannings  n  =  0.0238, 
as  found  by  Biery  (4),  was  verified  and  used  for  the  rough  boundaries. 
For  a  more  detailed  description  of  the  rough  boundary  pattern  and  the 
calibration  calculations  done,  see  reference  (20).  "Roughness  spacing 
in  Rigid  Open  Channels" • 

Geometries  I- a  and  I-b  were  tested  for  the  rough  boundary  condi- 
tion. The  results  indicated  that  a  similar  procedure  could  be  followed 
in  testing  of  other  geometries.  A  suitable  range  for  the  relative 
length,  —  ,  was  decided  upon  for  the  future  models.  The  channel 
width  ratio,  M,  had  to  be  chosen  so  that  submergence  could  be  obtained 
for  the  available  flow.  All  these  variables  were  selected  from  the  data 
obtained  from  the  testing  of  geometries  I-a  and  lib,  and  by  observa- 
tion of  the  submerged  flow  behavior. 

Besides  the  four  models  made  for  geometry  I-a  and  twelve  models 
made  for  geometry  I-b,  six  models  were  made  for  geometry  V-b,  three 
models  for  geometry  VI  and  four  models  for  geometry  VU.  Figure  25 
shows  the  models  that  were  tested  under  submerged  conditions.  The 
tables  (I  -  VII)  show  the  tests  performed  on  each  geometry  for  sub- 
mergence. 

Construction  of  Models 

The  two  dimensional  models  were  made  from  3/^"  thick  marine 
plywood  and  painted  to  keep  them  from  rotting.  Due  to  the  high 
pressure  exerted  by  the  backwater,  the  models  were  stiffened  by  mount- 
ing two  angle  irons  on  the  downstream  face  of  each  model  as  shown  in 
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Figure  12,  detail  at  test  section.  Holes  for  fastening  purposes 
were  made  along  the  vertical  sides  of  the  model. 

The  three  dimensional  models  were  made  from  3/8"  thick  exterior 
plywood.  These  models  were  braced  from  the  inside  by  2"  x  4M  studs. 
The  plywood  pieces  making  up  a  model  were  glued  together  to  prevent 
leakage  of  water  into  the  model.  The  three-dimensional  openings  were 
lined  with  sheet  metal.  A  special  DAP  Architectural  Grade  Caulking 
compound  that  could  adhere  readily  to  both  wood  and  sheet  metal  was 
used  as  seal  around  the  model  openings.  The  models  were  also  painted- 
to  prevent  the  plywood  from  rotting. 

Selection  of  Variables 
The  29  selected  models  were  to  be  tested  in  the  6k   foot  long 
flume,  first  with  smooth  and  then  with  rough  boundaries.  The  large 
number  of  variables  involved  made  it  necessary  to  establish  a  systema- 
tic procedure  for  selecting  the  quantities  allowed  to  vary  and  those 
held  constant  in  the  several  tests  so  that  the  possible  range  of  the 
variables  could  be  covered  uniformly  and  completely.  This  was  of 
concern  because  the  range  of  discharges  and  slopes  for  which  sub- 
mergence could  be  attained  was  very  limited  and  each  test  was  time 
consuming.  The  procedure  followed  to  solve  this  problem  was  based 
on  the  theoretical  analysis  and  on  the  testing  of  simple  geometries. 
The  theoretical  analysis  had  resulted  in  a  dimensionless  equation 
(24),  that  could  be  interpreted  to  prevent  the  grouping  of  the  data. 
The  testing  of  simple  geometries  also  gave  the  possible  range  of  the 
key  variables:  discharge,  slope  and  opening  ratio. 
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Preliminary  test  runs  had  given  information  about  the  channel 
opening  ratios,  M  =  b/B,  which  could  be  submerged  with  the  available 
flow  in  the  laboratory.  Testing  of  the  simple  geometry  I-a  with  the 
selected  channel  opening  ratios  gave  the  flow  rates  necessary  for 
each  model  to  have  submergence.  Now,  by  interpreting  the  dimension- 
less  equation  (2*0  the  ordinate s  of  the  plotting  points  on  the  graphs 
could  be  computed  for  the  several  flow  rates  of  the  respective  channel 
opening  ratios  chosen.  A  table  was  made  containing  these  ordinates 
for  five  discharges  '  and   the  opening  width  representing  the  five 
channel  widths  ratios  to  be  investigated.  The  dimensionless  ordinate 
axis  was  found  to  be: 


k*-' 


by  substituting  the  different  values  of  Q  (cfs)  and  b  (ft),  g  being 
constant,  the  following  table  was  obtained  for  the  ordinate  values 
for  the  five  channel  opening  ratios,  M  ," 


Q 

A'/* 

cfs 

M'    =  0.15 

M-    =  0. 

20 

M!    =  0.25 

M!   =  0. 

275 

M1    =0.30 

2.0 

0.7095 

0.353 

0.202 

0.159 

0.128 

2.5 

0.885 

0.441 

0.2525 

0.199 

0.160 

3.0 

1.062 

0.529 

0.303 

0.238 

0.192 

3.5 

- 

- 

0.354 

0.278 

0.224 

4.0 

- 

- 

0.405 

0.313 

0.256 
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By  choosing  flow  rates  in  the  test  runs  that  give  different  ordinate 
values  in  the  table,  grouping  of  data  in  the  graphs  can  be  prevented. 
As  an  example  it  can  be  seen  from  the  table  that,  to  give  a  complete 
uniform  coverage  for  plotting  of  a  curve  it  is  not  necessary  to  make 
more  than  about  k  test  runs  if  they  axe  carefully  selected.  When 
using  the  table,  it  should  be  kept  in  mind  that  it  would  not  be 
possible  to  make  all  21  test  runs  listed  in  the  table  for  all  geo- 
metries even  if  desired.  This  is  because  some  test  geometries  may 
be  overtopped  and  some  will  not  obtain  submergence,  depending  on 
the  value  chosen  for  the  Froude  number.  The  seven  tables  that  follow 
this  section  show  the  test  runs  made  for  the  different  geometries. 
The  conditions  selected  in  Table  I  through  VII  were  chosen  in  part 
from  the  table  developed  in  this  section.  Further  limitations  were 
imposed  as  to  number  of  test  runs  after  the  data  taken  for  the 
limiting  curves  were  calculated  and  plotted.  From  the  spacing  of 
these  curves,  the  decision  was  made  as  to  the  number  of  additional 
runs  desired.  If  the  limiting  curves  were  closely  spaced,  runs  were 
made  to  establish  how  the  other  curves  were  spaced  between  the 
limiting  ones. 
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Table  I  -  Smooth  Boundary  Tests 

n  =  0.0110 

Geometry  I-a 

Flow 

Froude  Number 

Channel  Width 

Ratio 

(b/B) 

Rate 

Fn 

• 

cfs. 

0.10    0.15   0.20 

0.20 

0.25 

0.275 

0.30 

2.0 

X 

X 

X 

X 

X 

2.0 

X 

X 

* 

* 

* 

2.0 

X 

X 

* 

* 

* 

2.5 

X 

0 

X 

X 

X 

2.5 

X 

X 

X 

X 

X 

2.5 

X 

X 

X 

X 

X 

3.0 

X 

0 

X 

X 

X 

3.0 

X 

0 

X 

X 

X 

3.0 

X 

0 

X 

X 

X 

3.5 

3.5 

X 

0 

X 

X 

X 

3*5 

X 

0 

X 

X 

X 

X  Regular  test  measurements  obtained 

0  The  flow  overtopped  the  model  and  Chanel  walls 

*  Submergence  was  not  possible  for  these  flow  conditions 
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Table 

1  II 

-  Smooth  Boundary  Test 

• 

n 

=  0. 

011C 

I 

Geometry  I-b 

Flow 
Rate 
cfs. 

i 
Froude  Number 
Fn 
.10  .15  .20 

Channel  Width  Ratio  b/B 
0.20         0.25 
Thickness  Factor  L/b 
.25  .50  .75  1.0  .25  .50  .75 

1.0 

.25 

0. 
.50 

30 
•  75 

1.00 

2.0 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

2.0 

Xi 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

2.0 

Xi 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

2.5 

X 

0 

0 

0 

0 

X 

X 

X 

X 

2.5 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

2.5 

X 

X 

X 

X 

X 

X 

X 

X 

3.0 

Xi 

0 

0 

0 

0 

3.0 

X 

0 

0 

0 

0 

X 

X 

X 

X 

X 

X 

X 

X 

3.0 

X 

0 

0 

0 

0 

X 

3.5 

Xi 

0 

0 

0 

0 

0 

0 

0 

0 

3.5 

X 

0 

0 

0 

0 

X 

X 

X 

X 

X 

X 

X 

X 

3.5 

X 

0 

0 

0 

0 

X 

X 

4.0 

Xi 

0 

c 

0 

0 

0 

0 

0 

0 

4.0 

X 

0 

0 

0 

0 

0 

0 

0 

0 

X 

X 

4.0 

xi 

0 

0 

0 

0 

X  Regular  test  measurements  obtained 

0  The  flow  overtopped  the  model  and  channel  walls 

*  Submergence  was  not  possible  for  these  flow  conditions 

X.  Experimental  condition  investigated 
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Table  III  -  Rough  Boundary  Tests 
n  =  0.0238 
Geometry  I-a 


Flow 
Rate 
cfs. 

Froude  Number 
Fn 
0.10   .015    0.20 

2.0 

X 

2.0 

X. 

1 

2.0 

X 

2.5 

2.5 

X 

2.5 

X 

3.5 

3.5 

X 

3.5 

X 

Channel  Width  Ratio  (b/B) 
0.20     0.25     0.275   0.30 


X  Regular  test  measurements  obtained 

0  The  flow  overtopped  the  model  and  channel  walls 

*  Sumbergence  was  not  possible  for  these  flow  conditions 

X.  Experimental  condition  investigated 
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Table  IV  -  Rough  Boundary  Tests 
n  =  0.0238 
Geometry  I-b 


Froude  Number  Channel  Width  Ratio  (b/B) 

Flow      Fn  0.20         0.25  0.30 

Rate   .10  .15  .20  Thickness  Factor  (L/b) 

cfs.  .25  .50  .75  1.0  .25  .50  .75  1.0  .25  .50  .75  1.0 


2.0 
2.0 
2.0 
2.5 
2.5 
2.5 
3.5 
3.5 
3.5 


X 


X. 

1 


xi 


X 

X 

X 

X 

* 

* 

0 

0 

0 

0 

X 

X 

X 

X 

X 

X 

X 

X 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

X 

X 

xi 

0 

0 

0 

c 

X  Regular  test  measurements  obtained 

0  The  flow  overtopped  the  model  and  channel  walls 

*  Submergence  was  not  possible  for  these  flow  conditions 

X.  Experimental  condition  investigated 
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Table 

V  . 

-  Rough  Boundary  Tests 

n  =  0. 

0238 

Geometry  V-b 

Flow 
Rate 
cfs. 

Froude  1 

Fn 

.10       .15 

dumber 
.20 

Thickness  Factor  L/b 

Channel  'Width  Ratio   i 

0.20 

Angle  of  Skew0    2 

15°     22|°     30°       15° 

=  0.50 

(b/B) 

0.30 

22i° 

30° 

2.0 

X 

X         X 

X           X 

X 

X 

2.0 

X 

X         X 

X           * 

* 

* 

2.0 

X 

X         X 

X            * 

* 

* 

3.0 

X 

X 

X           X 

X 

3.0 

X 

X 

X           X 

X 

3.0 

X 

X 

X            X 

X 

X  Regular  test  measurements  obtained 

0  The  flow  overtopped  the  model  and  channel  walls 

*  Submergence  was  not  possible  for  these  flow  conditions 
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Table  VI 

-  Rough  Boundary 

Tests 

n  = 

0.0238 

Geometry  VI 

Flow 

Froude 

Number 

Thickness 

Factor 

L/b  - 

-   0.50 

Rate 

Fn 

Channel  Width  Ratio 

b/B 

cfs. 

.10 

.15 

.20 

.15 

.20 

.25 

2.0 

X 

X 

X 

X 

2.0 

X 

X 

X 

* 

2.0 

X 

X 

X 

* 

2.5 

X 

X 

X 

X 

2.5 

X 

X 

X 

X 

2.5 

X 

X 

X 

* 

3.0 

X 

X 

X 

X 

3.0 

X 

X 

X 

X 

3.0 

X 

X 

X 

X 

X  Regular  test  measurements  obtained 

0  The  flow  overtopped  the  model  and  channel  walls 

*  Submergence  was  not  possible  for  these  flow  conditions 
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Table  VII  -  Rough  Boundary  Tests 

n  =  0.0238 

Geometry  VII 

Flow   Froude  Number  Channel  Width  Ratio  b/B 

Rate       Fn  0.35  0.50 

cfs.  .1  .2  .3  A   .5  (3=   .0  /5=  .3  /?=  .5  j3  =   -0  /?  =  .3  /3  =  -5 


*  *  X 

*  *  X 

*  *  X 

*  *  X 

*  *  X 

*  * 


1.0 

X 

* 

2.0 

X 

* 

X 

2.0 

X 

* 

X 

2.0 

X 

* 

X 

2.0 

X 

* 

2.0 

X 

* 

2.5 

X 

X 

2.5 

X 

* 

X 

2.75 

X 

* 

X 

3.0 

X 

3.0 

X 

X 

3.85 

X 

3.85 

X 

3.85 

X 

X 

3.85 

X 

X 

3.85 

X 

X 

0  X 

0  *  X 

0  *  *      X 

0  *  *       X 

0  *  *      X 


X  Regular  test  measurements  obtained 

0  The  flow  overtopped  the  model  and  channel  walls 

*  Sumbergence  was  not  possible  for  these  flow  conditions 
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CHAPTER  VI 
TESTING  PROCEDURE  AND  DATA  PROCESSING 

Test  Setups  for  Submerged  Flow  Conditions 

In  the  previous  chapter  the  geometries  to  be  tested  were  described 
and  the  selection  of  variables  were  made.  To  set  up  a  test-condition 
where  uniform  flow  exists  with  a  pre-determined  discharge  and  Froude 
number  in  the  unobstructed  channel  can  be  a  rather  time-consuming 
process  when  it  is  based  upon  a  trial  and  error  method.  P.  F.  Biery 
(4)  had  developed  a  systematic  procedure  to  obtain  the  uniform  flow 
in  the  flume  after  a  Froude  number  and  a  flow  rate  had  been  selected. 
A  similar  procedure  was  adopted  for  the  present  investigation. 

Nomographic  Charts 

A  nomographic  chart  was  prepared  which  relates  the  Froude  number 
F  ,  the  normal  depth  Y  ,  the  discharge  Q,   the  tailgate  elevation  Y.  , 
the  channel  slope  S,  and  the  slope  control  reading.  The  relationships 
used  are: 


Fn  "  BY„  r=-  (25) 

Manning's  formula 

v  =  1^86    r2/3s1/2  (26) 
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Francis  weir  formula 


Q  =  K  B 


(H  +  H  )3/2  -  H  3y/2 
v     v'      V 


(27) 


The  ctxrves  relating  these  variables  are  shown  in  Figures  26 
and  27  for  smooth  and  rough  boundaries  respectively.  The  calcula- 
tions were  made  using  the  1620  IBM  computer.  The  program  is  given 
together  with  the  derivation  of  the  relationships  in  Appendix  C. 

To  use  the  monographic  chart  one  enters  the  left  side  with  a 
selected  Froude  number  and  discharge  and  read  the  normal  depth.  Now, 
with  Q  and  Y  fixed,  the  slope  is  found.  With  Q  and  slope  fixed, 
the  tailgate  setting  Y.  is  read.  The  slope-revolution  calibration 
curve  is  projected  onto  the  nomograph  and  the  setting  of  the  revolution 
counter  can  be  read  directly  from  the  charts.  Arrows  drawn  on  the 
charts  demonstrate  the  path  to  follow  for  the  selection  of  a  particular 
test  condition.  The  use  of  these  charts  helped  to  set  the  normal 
depth  flow  within  30  minutes  even  when  corrections  had  to  be  made. 
The  normal  depth  was  always  checked  along  the  centerline  and  minor 
adjustments  made  of  the  tailgate  to  ensure  uniform  flow  in  the  flume. 

After  the  test  conditions  were  determined  and  set  by  means  of 
the  chart  of  fig.  26, ,  :he  angle  irons  were  fastened  at  the  test  section 
and  the  model  was  placed  behind  them  to  the  correct  position  and  bolted 
to  the  angle  irons.  The  backwater  rose  behind  the  model  until  equili- 
brium was  obtained  between  the  upstream  and  downstream  flow.  The 
time  required  to  obtain  equilibrium  varied  between  10  to  30  minutes 
depending  on  the  opening  of  the  model,  the  Froude  number  and  the  rate 
of  flow  used. 
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FIGURE  2  6  TESTS    SELECTION    CURVE  —  LARGE     FLUME  —  SMOOTH    BOUNDARES 
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FIGURE  27   TESTS  SELECTION   CURVE 


LARGE   FLUME  —  ROUGH   BOUNDARIES 
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Figure  28  shows  an  upstream  and  a  downstream  view  of  a  typical 
test  run  after  the  backwater  has  come  to  equilibrium.  At  this  stage, 
surface  elevation,  streamlines,  velocities  and  flow  patterns  were 
observed  and  recorded.  All  depths  were  measured  relative  to  the 
bottom.  The  electronic  point  gage  was  used  to  measure  the  backwater. 
Dye  was  injected  by  means  of  a  syringe  having  an  18"  long  needle,' 
in  order  to  study  the  streamline  paths  and  flow  behavior  as  presented 
in  Figures  3  and  38.  The  Prandtl  tube  connected  to  the  pressure 
diaphragm  was  used  for  all  velocity  readings  (see  section  on  measur- 
ing equipment.  Chapter  IV). 

Data  Processing 
The  data  obtained  from  the  test  runs  were  processed  on  the 
134  1620,  IBM  7090  and  LGP  30  computers.  By  using  the  computers, 
valuable  time  was  saved,  and  more  reliable  answers  were  obtained 
than  could  be  expected  from  manual  calculations.  An  outline  of  the 
computer  programs  is  presented  in  Appendix  C. 
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FIG.    28  UPSTREAM       AND       DOWNSTREAM       VIEWS      OF 

TYPICAL      TEST     RUN 
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CHAPTER  VII 
RESULTS  AND  THEIR  ANALYSIS 

In  this  chapter  the  results  of  the  processed  data  are  presented 
and  analyzed.  Due  to  the  large  amount  of  data  taken,  the  test  runs 
will  not  be  treated  separately  but  in  groups  according  to  the  geo- 
metries.  The  results  from  each  group  of  tests  are  presented  in 
dimensionless  graphs  and  an  explanation  is  given  as  to  each  graphs 
significance. 

The  terms  free  outflow,  fully  submerged  and  partly  submerged 
outflow  are  frequently  used  in  this  chapter  and  are  defined  as  follows. 
Free  outflow  exists  when  the  jet  discharging  from  the  downstream  face 
of  the  bridge  opening  is  surrounded,  except  for  the  channel  bottom, 
only  by  air  and  a  hydraulic  jump  occurs  further  downstream.  Fully 
submerged  outflow  exists  when  the  jet  discharges  into  a  stream  that 
has  a  depth  equal  to  or  greater  than  the  downstream  height  of  the 
bridge  opening.  A  partly  submerged  outflow  is  the  case  when  the  jet 
discharges  into  a  stream  that  has  a  depth  less  than  the  height  of 
the  bridge  opening  and  no  downstream  jump  occurs. 

The  coefficients  of  contraction,  of  discharge  and  of  velocity 
are  presented  for  geometry  I-a  with  free  jet  outflow.  For  the  same 
flow  conditions,  the  velocity  distribution  over  the  cross  section  at 
vena  contracta  is  presented  in  form  of  isovelocity  curves  covering 
the  cross  section.  A  velocity  comparison  at  different  elevations 
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in  the  vena  contracta  is  also  given  some  special  attention. 

The  backwater  depths  obtained  for  the  different  geometries  are 
presented  first  in  terms  of  the  discharge  and  with  the  Froude  number 
as  a  parameter,  and  second  as  a  function  of  the  ratio  of  the  Froude 
number  to  the  channel  opening  ratio. 

The  presentation  of  the  data  can  be  broken  into  two  major 
categories  -  the  smooth  boundary  (n  =  0.0110)  and  the  rough  boundary 
(n  =  0.0238).  Comparisons  are  made  between  findings  from  the  two 
boundary  conditions.  Pictures  of  typical  directions  of  the  dis- 
charge jet  for  different  flow  conditions  are  presented  together  with 
a  discussion  of  other  observations  and  analysis  made. 

Geometry  I-a  -  Smooth  Boundaries 

Coefficients  of  Contraction,  of  Discharge  and  of 

Velocity  for  Free  Jet  Outflow 

In  order  to  investigate  the  coefficients  of  velocity,  contrac- 
tion and  discharge  for  submerged  bridge  inlets  6?  test  runs  were  made 
(see  Table  I  -  VH).  The  coefficients  are  defined  as  follows: 
Coefficient  of  contraction,  C  ,  is  the  ratio  of  jet  area  at  vena 
contracta  A  to  the  area  of  the  orifice  A_. 


0 


c   A 
o 


Coefficient  of  discharge,  C,t  is  the  ratio  of  the  actual  dis- 
charge Q  ,  obtained  from  the  Venturineter  reading,  to  the  theoretical 
Q. ,  as  calculated  by  formula  (10) 
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Coefficient  of  velocity  is  the  ratio  of  actual  velocity  V  to  the 

theoretical  velocity  V, 

V 
c  ■  -& 

V  \ 

It  is  customary  to  combine  the  three  coefficients  as  follows: 

From  continuity  equation  the  following  can  be  written: 

Q    A  V 
a    v  a 


Qt   A.  V. 


o  t 
from  which  it  follows  that: 


d    c  v 


To  find  the  coefficient  of  contraction  the  vena  contracta  had 
to  be  located  and  the  cross-sectional  area  measured.  This  was  done 
by  means  of  the  electronic  pointgage. 

A  sufficient  number  of  readings  were  taken  at  the  vena  contracta 
to  plot  the  cross  sectional  area  on  graph  paper.  By  use  of  a  plani- 
meter  the  area  A  was  obtained  and  the  coefficient  of  contraction 
could  be  calculated,  as  the  area  of  the  model  opening  was  known. 
The  coefficient  of  contraction  was  plotted  versus  the  ratio  b/Y. 
of  the  model  opening  width  to  the  backwater  depth  using  the  channel 
width  ratio  b/B,  as  parameter  (see  Figure  29). 

To  obtain  values  for  the  calculation  of  the  coefficient  of  dis- 
charge C,,  a  theoretical  discharge  equation  was  developed.  The 
development  of  this  equation  (10)  was  presented  in  the  theoretical 
analysis  in  Chapter  m.  A  detailed  presentation  of  the  numerical 
calculations  involved  is  given  in  Appendix  B.  The  actual  rate  of 
flow  through  the  constriction  was  obtained  from  the  manometer  readings 
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and  the  venturi  meter  calibration  curves.  Curves  were  plotted  for 
the  coefficient  of  discharge  C.  versus  the  ratio  b/Y..  and  the  channel 
width  ratio  was  again  used  as  parameter  (see  Figure  29).  The  co- 
efficient of  velocity  was  then  found  from  the  relationship  between 
the  three  coefficients,  C.  =  C  C  and  plotted  versus  b/Y,  in 
Figure  29. 

The  C.  curves  seem  to  converge  toward  the  value  0.61  as  the  -k- 
ratio  tends  to  zero.  The  curves  are  closely  spaced  straight  lines 
and  the  C,  values  decrease  with  increasing  channel  width  ratios. 

The  C  curves  also  tend  towards  the  0.61  value  as  the  -st*  ratio 
c  Y1 

decreases.  Values  of  C  decrease  for  increasing  channel  width  ratios. 

For  increasing  C  values,  within  the  test  range,  the  C  curves  tend 

to  increase  toward  the  value  0.71  and  then  stabilize  at  this  value. 

The  C  curves  have  the  same  relative  spacing  as  the  C  and  C,  curves 

as  to  the  —  ratio  but  tend  toward  the  value  1.0  as  a  limiting  value 

when  the  7  ratio  decreases  toward  zero. 
x1 
These  findings  check  well  with  the  typical  value  0.611  found 

21 
from  the  theory  of  free  streamlines.  Streeter     shows  this  for  a 

two  dimensional  ideal  orifice  where  he  applies  complex  variables  to 
the  "Schwarz  -  Christoffel  Transformation"  and  gets  the  ideal  dis- 
charge coefficient 


n  2b  IT  , 

77 

22 
Nearly  a  century  ago  Kirchhoff     showed  that  the  contraction 

coefficient  for  irrotational  efflux  from  a  plane  orifice,  under 

idealized  conditions,  would  have  the  magniture   ,  ?   =  0.611. 
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1  3 
Von  Mises     determined  tables  of  the  C  coefficient  for  several 

c 

boundary  characteristics  which  show  the  same  trend  as  those  found 

in  this  investigation.  The  test  data  are  recorded  in  test  runs 

1  -  67  shown  in  Tables  (VIII  -  XV).  The  calculated  values  to  find 

C  .  C.,  and  C  are  found  in  Table  (XXX). 
c   d      v 

Velocity  Distribution  at  Vena  Contracta  for  Free  Jet  Outflow 
Extensive  velocity  measurements  were  taken  in  the  jet  downstream 
of  constriotion  at  the  vena  contracta  in  order  to  observe  the  velocity 
distribution.  The  Prandtl  tube  together  with  the  pressure  transducer 
were  used.  The  readings  were  taken  in  a  grid  system  to  obtain  full 
coverage  over  the  cross  section.  The  borderline  of  the  area  at  the 
cross-section  was  also  measured  and  plotted  on  graph  paper  together 
with  the  velocity  grid  values.  From  this  the  isovelocity  curves 
were  drawn.  Figure  30  shows  one  of  the  typical  isovelocity  curves 
drawn  from  such  grid  measurements.  Figure  ^\    shows  the  result  obtained 
by  following  the  individual  isovelocity  curves  in  the  water  cross- 
section  itself.  This  was  done  by  moving  the  Prandtl  tube  through 
the  cross-section  and  mapping  the  path  for  the  individual  isovelocity 
curves  obtained  by  means  of  the  horizontal  and  vertical  traversing 
devices.  The  data  of  Figures  30  and  31  correspond  to  discharges 
3.75  cfs  and  2.62  cfs,  normal  depth  Froude  number  0.1*4-2  and  0.103  and 
opening  ratios  0.30  and  0.30  respectively. 

In  Figure  32  the  values  of  the  velocities  are  plotted  at  different 
depths  and  a  mean  velocity  curve  is  drawn.  The  data  for  Figure  32 
are  the  same  as  those  for  Figure  31 •  Figure  32  shows  the  velocity 
variation  over  the  opening  width  at  different  depth.  A  similar 
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study  was  made  for  other  flow  conditions  as  recorded  in  Tables  XXVIU 
and  XXIX.  All  the  recorded  values  when  plotted,  show  the  same 
characteristic  velocity  distribution  as  shown  in  Figures  30,  31. 
and  32. 

The  Generalized  Backwater  Equation  - 
Outflow  Fully  and  Partly  Submerged 
From  the  theoretical  analysis  done  in  Chapter  III  it  is  expected 
that  the  general  relation  between  the  backwater  and  normal  depth  is 
of  the  form 


Y      C 
n 


<  #   )  I  (36) 


where  C  is  a  coefficient  that  would  take  into  account  the  effects 
of  the  approach  velocity,  discharge  coefficient,  non-hydrostatic 
pressure  distribution,  non-uniform  velocity  distribution  and  other 
empirically  determined  factors.  In  a  logarithmic  plot  of  Y>/Y 
versus  (Fn/M»)  ,  the  slope  of  the  line  is  / .  A  total  of  33  tests 
were  made  on  geometry  I-a  and  processed  on  the  IBM  1620  computer  to 
obtain  plotting  values  of  (Y, /Y  -  1)  versus  (F  /K1)  .  The  value 
(Y.  /Y  -  1 )  was  used  in  order  to  expand  the  scale  of  the  graph. 
The  graph  is  shown  on  Figure  33»  From  this  graph  the  values  y  and  C 
were  round  to  be  0.90  and  1.18  respectively,  and  equation  (36)  can  be 
written  for  the  submerged  conditions  of  geometry  I-a  in  the  following 
way: 

0.90 
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or 

J>    =     1  +1J8  (  $   )1'80  •        (18) 

n 

The  test  data  taken  to  obtain  this  generalized  backwater  equation  (18) 

are  recorded  in  test  runs  102-1 3^*-  which  are  located  in  Tables  XVI 

and  XVH.  The  calculated  values  from  the  data  are  presented  in 

Table  XXX. 

A  general  backwater  equation  was  not  obtainable  for  the  free 

jet  outflow.  This  is  because  the  Froude  numbers  for  test  runs  1  -  67 

are  all  in  the  range  of  0.550  -  0.587  which  re  stilts  in  one  group  of 

numbers  and  a  line  could  not  be  defined.  The  calculation  of  the  Froude 

numbers  on  the  computer  are  shown  in  Appendix  C  and  the  results  are 

listed  in  Table  XXX. 

Discharge  Coefficient  for  Free.  Submerged  and 
Partly  Submerged  Outflow 
A  graph  was  made  of  the  discharge  coefficient,  C,,  versus  the 
channel  opening  ratio,  M't  for  the  submerged  smooth  bridge  opening, 
Geonetry  I-a  (Figure  3^)«  The  curves  obtained  are  for  the  submerged 
discharge  jet  except  for  the  two  upper  most  curves  which  are  for 
partly  submerged  jets.  These  curves  are  similar  to  those  estimated 
by  P.  F.  Biery  (k)   for  the  unsubmerged  case  and  for  channel  opening 
ratios  of  0.29  and  up.  The  present  values  for  the  submerged  case 
were  obtained  for  the  range  of  M*  between  0.10  and  0.35«  Figure  J>h 
shows  the  C,  versus  M'  curves  for  submerged  and  partly  submerged  flow 
as  well  as  Biery* s  discharge  coefficient  curves  for  the  unsubmerged 
case  which  have  been  superposed  for  the  purpose  of  comparison.  The 
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test  data  taken  to  obtain  the  curves  in  Figure  3^  are  recorded  in 
test  runs  102  to  134  which  are  located  in  Tables  XVI  and  XVII.  The 
calculated  values  are  presented  in  Table  XXX. 

The  values  of  the  discharge  coefficient  are  summarized  in  Fig- 
ure  35»  For  the  submerged  jet,  two  sets  of  curves  are  given,  one 
uses  the  Froude  number  as  a  parameter,  the  other  uses  the  ratio 
Y  /  r  of  the  normal  depth  to  the  arch  radius.  The  free  jet  discharge 
coefficient  curves  previously  presented  in  Figure  29  have  been  re- 
produced here  for  completeness.  The  calculations  are  shown  in 
Table  XXX.  . 

Geometry  I-b  -  Smooth  Boundaries 
As  can  be  seen  from  the  definition  of  test  geometries,  the 
geometry  I-a  is  a  special  case  of  geometry  I-b.  By  setting  the 
variable  r-  equal  to  zero  in  the  definition  of  geometry  I-b,  the 
geometry  I-a  is  obtained.  Geometry  I-a  therefore  represents  the 
limiting  case  of  geometry  I-b,  and  can  be  treated  together  with  I-b. 
In  the  theoretical  analysis  made  in  Chapter  III  the  dimension- 
less  relation 

gV2b5/2  =  f  (  T  } 

was  developed  for  the  orifice  type  flow  through  a  submerged  arch 
constriction.  A  total  of  92  test  runs  were  conducted  for  submerged 
conditions  in  the  smooth  channel.  Three  sets  of  dimensionless  curves 
were  plotted.  To  save  space  and  make  it  easier  to  compare  the  curves, 
they  were  all  plotted  in  one  figure.  The  r  ratio  is  used  as  para- 
meter, and  the  Froude  number  is  constant  for  each  set  of  curves. 
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To  be  able  to  show  all  three  sets  of  curves  clearly  on  one  graph, 
the  vertical  coordinate  axis  was  broken  into  three  parts  each  start- 
ing from  zero.  The  curves  are  shown  in  Figure  J6,     For  higher 
Froude  numbers  the  curves  get  closer  spaced.  This  means  that  the 
L/b  ratio  becomes  insignificant  for  all  practical  purposes  when  the 
Froude  number  reaches  a  value  of  about  0.20  for  the  smooth  (n  =  0.0110) 
boundary  case  of  geometry  I-b.  The  data  from  the  92  test  runs  are 
listed  in  run  number  101  -  192  and  recorded  in  Tables  XVI  -  XXI. 
The  calculated  values  of  the  plotted  points  are  shown  in  Tables  XXX 
and  XXXI. 

Geometries  I-a  and  I-b  -  Rough  Boundaries 
After  the  roughness  was  installed  and  calibrated  for  Mannings 
n,  test  selection  curves  were  computed  and  plotted  as  shown  in 
Figure  27.  An  average  Mannings  n  of  0.0238  was  found  and  used  for 
the  rough  boundary  condition. 

Geometries  I-a  and  I-b  were  again  treated  together.  A  total  of 
38  test  runs  were  made  in  order  to  plot  the  dimensionless  curves 
shown  in  Figure  37  •  These  three  sets  of  curves  were  also  based  on 
the  dimensionless  equation  (24)  developed  in  Chapter  m.  Due  to  the 
previous  experience  with  Geometries  I-a  and  I-b  in  the  smooth  boundary 
test  runs,  a  better  test  selection  was  possible  for  the  rough  boundary 
tests  with  respect  to  grouping  of  the  points  in  the  graphs.  The 
data  are  recorded  in  test  runs  193  -  231  and  can  be  found  in  Tables 
XXI  -  XXIII.  The  calculated  values  from  the  data  are  listed  in 
Table  XXXI. 
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The  investigation  shows  that  the  width  of  the  constriction  in 
the  direction  of  flow,  L,  influences  the  backwater  to  a  certain 
degree  when  the  Froude  number  is  low.  The  wider  the  bridge  in  the 
direction  of  flow  the  lower  is  the  backwater  depth  obtained.  For  a 
Froude  number  of  0.15  the  backwater  ratio  for  an  L/b  ratio  of  1.0 
was  as  much  as  30$  lower  than  for  a  L/b  ratio  of  0.00.  For  the  more 
practical  cases  where  the  Y./b  ratio  is  in  the  order  of  about  0.50 
to  1.00,  a  lowering  of  about  10  -  20$  in  the  backwater  ratio  is 
usual  when  the  L/b  ratio  is  increased  from  0.00  to  1.00.  Flow  through 
three  dimensional  arch  bridges  is  very  similar  in  behavior  to  the  flow 
through  short  culverts.  For  short  culverts,  the  flow  characteristics 
are  relatively  independent  of  the  slope  and  the  factors  that  involve 
the  length  are  comparatively  unimportant.  The  control  section  for 
square  edge  openings,  such  as  investigated  here,  is  at  the  inlet. 

Separation  of  flow  around  the  opening,  at  the  inlet  section, 
was  frequently  observed.  The  flow  contracted  towards  vena  contracta 
section  in  the  first  part  of  the  barrel  expanded  in  a  spiral  motion 
as  shown  in  Figure  38  and  ended  in  a  slug  flow  at  the  submerged  out- 
let. Figure  39  shows  the  case  were  slug  flow  occurs  at  the  exit  of 
the  barrel.  Figure  40  presents  the  free  jet  case  for  geometry  I-a 
with  rough  boundaries  and  the  same  flow. 

Vortex  action  was  also  often  observed.  This  air-vortex  occurs 
from  the  surface  right  above  the  crest  of  the  inlet  and  introduces 
air  through  the  opening.  This  decreases  the  water  area  of  the  opening 
and  hence  the  flow  through  the  constriction,  and  as  a  result  the 
backwater  increases. 


94 


95 


Fig  38a   Spiral    Motion    in    Barrel    Section    Downstream    of 
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Fig  38b     Typical     Flow     Condition    Through     Constriction 
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FIG.      39  SLUG        FLOW       AT       BARREL      EXIT 


FIG      40  FREE  DISCHARGE         JET 
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An  important  effect  upon  the  separation  phenomena,  that  controls 
the  change  of  flow  from  a  sluice  type  to  full  conduit  flow,  is  the 
existence  of  a  negative  pressure  in  the  barrel.  This  is  caused  by 
the  relative  magnitudes  of  air  inflow  from  the  outlet  and  from  the 
vortex  action  over  the  inlet.   The  outflow  of  air  due  to  entrairnent 
at  the  air-water  surface  also  has  an  effect  upon  the  separation 
phenomena,  the  more  air  entrained  the  less  separation  will  occur  and. 
the  barrel  will  tend  towards  full  conduit  flow. 

Ccr.r-arison  Between  Smooth  and  Rou^h  Boundary  Tests 
A  comparison  between  the  limiting  curves  of  Figures  3^  and  37  for 
the  smooth  and  rough  boundary  conditions  respectively  are  made  in  Fig- 
ure 41 .  These  curves  show  that  a  rough  boundary  condition  is  of  advantage 
for  higher  Froude  numbers  (Fn>C.15),  as  co-old  be  exoected  from  observing 
the  flow  behavior.   For  higher  Froude  numbers  the  approach  velocities 
will  increase  and  consequently  give  a  higher  turbulence  level  due  to  the 
natural  roughness.   This  will  improve  the  ability  of  the  barrel  to  pro- 
duce full  conduit  flow,  because  it  increases  the  air  entrainment  in  the 
barrel  at  the  air-water  surface  and  hence  reduces  the  separation  at.  the 
inlet.   To  further  help  the  air  entrainment,  a  roughening  of  the  inlet 
surface  would  be  desirable.  It  appears  that  by  increasing  the  upstream 
approach  flow  turbulence  or  introducing  a  higher  turbulence  in  the  inlet 
section,  a  more  desirable  flow  capacity  is  obtained  through  the  con- 
striction and  hence  lower  backwater  depth.  In  other  words,  the  ex- 
periments show  that  for  structures  flowing  full  such  as  culverts  the 
aim  should  not  be  to  design  approach  channels  of  highly  engineered 
smooth  boundaries  and  perhaps  percise  alignment ,  but  rather  to  con- 


98 


o 
o 

e 

CO 


E 
o 
<r> 
CD 


*•>      a 

c 

W        O 

a>     CD 

> 


>=^ 


3 


v> 

c 
g 
to 
c 
<u 

E 
b 


o 
a. 

E 
o 
o 


o 


99 


centrate  on  high  approach  velocities  and  agitating  the  turbulence 
in  the  flow.  This  way  more  air  entrainment  is  obtained  and  hence 
more  security  of  getting  a  full  conduit  flow  through  the  constriction. 
The  separation  effect  can  also  be  greatly  reduced  by  use  of 

rounded  inlet  sections,  tapered  transition  sections  or  equivalent 

23 
inlet  structures  as  shown  by  J.  L.  French  -' ."  .  He  also  conducted 

experiments  indicating  that  full  conduit  flow  in  corrugated  barrels 

occurred  at  appreciably  lower  pool  levels  than  in  smooth  wall  barrels 

with  the  same  inlet  structure.  These  observations  made  by  French 

shows  the  same  trend  as  those  made  by  the  author  in  the  present  study 

of  the  flow  behavior  for  submerged  arch  bridge  constrictions. 

To  improve  the  flow  through  submerged  constrictions  or  culverts  the 

ability  of  the  constriction  to  produce  full  flow  must  be  perfected. 

From  the  above  discussion  of  the  observations  it  can  be  summarized 

that  the  following  characteristics  are  of  essence. 

1 )  Reduction  of  separation  effects  through  the  inlet. 

2)  Changes  in  the  geometry  of  the  approach  channel. 

3)  Increase  in  air  entrainment  at  the  inlet  and  in  the  barrel. 
*0  Reduction  in  vortex  formation  above  the  inlet. 

5)  Reduction  in  air  inflow  into  the  barrel  from  the 
constriction  outlet. 
The  expected  solutions  to  some  of  these  essential  characteristics 
were  discussed  above.  It  was  not  the  scope  of  this  study  to  investi- 
gate these  suggested  improvements,  but  further  experimental  work 
could  with  profit  proceed  in  this  direction. 
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Skewed  Bridges  -  Geometry  V-b 
Rough  Boundaries  n  =  0.0238 

Figures  42  a  and  b  show  the  flow  at  the  outlet  of  a  left  and 
right  skewed  bridge.  The  discharge  jet  shoots  to  the  side  of  the 
flume  in  a  direction  approximately  perpendicular  to  the  face  of  the  , 
bridge. 

The  test  results  are  presented  in  the  dimensibnless  performance 
curves  in  Figure  43,  based  on  the  theoretical  equation  (24)  developed 
in  Chapter  III.  Three  sets  of  curves  are  plotted  using  the  Froude 
number  and  the  angle  of  skewness,  #2,  as  parameters.  The  width 
ratio,  L/b,  is  kept  constant  at  the  value  L/b  =  0.50»  The  performance 
curves  show  that  a  variation  in  the  angle  of  skew  between  0?  =  1  5 
and  $2  =   30  is  a  relative  unimportant  factor  compared  to  the  effect 
of  a  change  in  the  Froude  number.  To  be  noted  is  also  the  better 
performance  of  these  skewed  bridge  constrictions  for  low  Froude 
numbers,  (F  =0.10)  compared  to  the  <j>~  =   0.0  constrictions  of 
Geometry  I-b.  The  explanation  may  be  that  the  non-symmetry  of  the 
skewed  bridge  constriction  generates  more  turbulence  in  the  flow  at 
the  inlet  and  hence  a  flow  condition  closer  to  conduit  flow,  at  a 
lower  Froude  number.  As  the  Froude  number  increases  and  both 
Geometries,  I-b  and  V-b,  agitate  about  the  same  amount  of  turbulence 
in  the  flow,  the  difference  in  the  performance  curves  decreases. 
The  part  of  the  curves  having  a  Y..  /b  ratio  less  than  about  0.5  do 
not  pertain  to  the  submerged  inlet  cases,  but  are  merely  an  expected 
trend  of  the  curves  for  the  free  surface  flow.  A  total  of  24  test 
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FIG      42a  FLOW      THROUGH       LEFT        SKEWED     BRIDGE 


FIG      42b         FLOW       THROUGH        RIGHT        SKEWED      BRIDGE 
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runs  were  made  on  the  skewed  bridges.  They  are  recorded  in  Tables 
XXIII  and  XXIV  as  test  run  number  232  to  255.  The  calculated  data 
to  plot  the  curves  are  listed  in  Table  XXXI. 

Geometry  VI  -  Two  Span  Semicircular  Arch'  Bridge  Constriction 
Rough  Boundaries  n  =  0.0238 

This  Geometry  represents  about  11$  of  all  existing  arch  bridges 
in  the  state  of  Indiana^  The  great  majority  of  these  bridge  con- 
strictions have  a  thickness  factor,  L/b,  of  about  0.50,  which  was 
chosen  for  the  experimental  models  of  Geometry  VI. 

The  test  results  are  presented  in  the  dimensionless  graphs  shown 
in  Figure  kb.     To  be  noted  is  that  the  opening  width  of  the  models, 
b,  is  the  width  of  one  of  the  openings.  The  performance  of  the 
modes  in  Geometry  VI  is  approximately  the  same  as  for  Geometry  I-b 
of  same  thickness  factor.. 

A  total  of  27  test  runs  were  made  and  the  data  taken  are  recorded 
in  test  runs  no.  (256  -  282)  as  shown  in  Tables  XXIV  -  XXVI.  The 
calculated  values  are  presented  in  Table  XXXI. 

Geometry  VII  -  Segment  Arch  Bridge 
Rough  Boundaries  n  =  0.0238 
With  reference  to  Figures  4  and  5f  it  is  readily  seen  that  the 
dimensionless  relationship  of  the  discharge  as  a  function  of  the 
backwater  depth  and  the  width  of  the  constriction  holds  true  also 
for  the  segment  arch  openings.  The  curves  shown  in  Figure  ^5  are 
based  on  this  relationship  as  expressed  in  Equation  (24),  derived  in 
Chapter  III.  A  total  of  22  test  runs  were  done  and  recorded  in 
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Tables  XXVI  -  XXVII.  The  calculated  values  are  found  in  Table   XXXI. 

Figure  45  shows  that  the  spacing  between  the  curves  become 
very  small  as  the  Froude  number  increases  beyond  0.3.  This  means 
that  the  Froude  number  becomes  insignificant  as  a  parameter  when  the 
degree  of  turbulence  increases  beyond  a  certain  value  (  F  >  0.3). 

Head  Loss  Considerations 
In  the  theoretical  analysis  in  Chapter  III  an  expression  for  the 
additional  head  loss  coefficient  due  to  the  bridge  constriction  was 
found,  Equation  (35). 


V 

V02/2g 


(  to   )2   ,  ^0  }2 
V  A  ;    K   A.  ; 

y      1 


(35) 


This  equation,  which  holds  for  submerged  bridge  constrictions  is  similar 

2k 
to  the  one  derived  by  Izzard     for  free  surface  flow  through  a 

bridge  constriction,  and  K  varies  primarily  with 

a)  The  channel  opening  ratio  M' 

b)  Type  and  shape  of  bridge  constriction  such  as  inlet 
shape,  wing  wall  angle  $. ,  skewness  angle  $_  and 
thickness  factor,  L/b 

c)  Number  of  bridge  spans,  N 

d)  Eccentricity  of  bridge  opening,  e 

e)  Froude  number  F 

Other  variables  of  less  effect  may  however  have  a  cumulative 
effect  upon  the  head  loss  coefficient  depending  on  the  problem  at 
hand.  The  head  loss  coefficient  K  was  computed  on  the  IEH  7090  from 
Equation  (35).  The  corresponding  value  of  the  channel  opening  ratio 
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M*  was  also  computed,  and  by  a  least  mean  square  routine  the  slope 
of  the  line  of  the  plot  K  versus  M*  was  found  in  arithmetic  scale. 
A  significant  value  such  as  L/b,  02  and  N  was  used  as  parameters 
both  for  smooth  and  rough  boundary  tests.  The  plot  of  these  curves 
are  shown  in  Figures  46  to  S5» 

For  the  smooth  boundaries  a  sufficient  number  of  data  was  avail- 
able to  use  the  least  square  method  on  the  computer.  In  the  case  of 
rough  boundaries,  however,  the  data  available  for  each  parameter 
were  not  always  sufficient  and  some  of  the  lines  had  to  be  fitted 
by  eye. 

The  data  used  for  the  calculations  are  listed  in  run?  101  -  306 
presented  in  Tables  XVI  -  XXVII.  The  calculated  values  from  the 
computer  are  listed  in  tables  XXX  and  XXXI. 
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FIG.   47         HEAD     LOSS       COEFFICIENT     FOR     GEOMETRY 
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L 


Ib      SMOOTH       BOUNDARIES, 


0.50 


2.0 


0.8 


o 

- 

K  =    2.031    -   3.858  M1 

© 

o 
0s 

0     o 
D 

O   < 

\o 

i 

D   X 

0         > 

0 
FIG.    49 


0.1 


M' 


0.2 


0.3 


HEAD     LOSS     COEFFICIENT      FOR      GEOMETRY 


Ib        SMOOTH      BOUNDARIES,     -j=-    =    0.75 


112 


tL.d. 

n 

9 

9 

V 

K  =  1.969  -   4.695  M1 

1.8 

o 

V 

o 

? 

Vj 

1.4 

\       C\ 

D 

9  :  Points      Not 
Considered , 

\    \J 

o     \ 

Outlet      Not 

o 

Sufficiently 
Submerged 

1.0 

o 

3  \ 
\    ( 

D 

o 

\  ° 

n  6 

0.1 


0.2 


0.3 


M' 


FIG.    50  HEAD     LOSS     COEFFICIENT     FOR     GEOMETRY 


Ib  t    SMOOTH      BOUNDARIES, 


=    1.00 


\i 

—  =  0.00 

^.25 

\0.50 
\0.75 

\l.00 

0 
FIG.    5  1 


0.1 


M' 


0.2 


0.3 


SUMMARY      OF     HEAD      LOSS       COEFFICIENT 
CURVES     FOR     GEOMETRIES      Ia ,   a     lb 
SMOOTH      BOUNDARIES 


114 


FIG.    52        HEAD     LOSS     COEFFICIENT      CURVES       FOR      GEO' 
METRIES     Ia  ,   8    lb,    ROUGH     BOUNDARIES 
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General  Backwater  Equations  for  Geometries  I-a. 
I-b.  V-b.  VI  and  VII  -  Smooth  and  Rough  Boundaries 

The  theoretical  development  of  these  equations  was  done  in 
Chapter  III,  and  earlier  in  this  chapter  the  generalized  backwater 
equation 


t-< 


F 
V  M'  '. 


y 


(36) 


was  plotted  for  the  case  of  fully  submerged  discharge  jet  in  the  form 
(see  Figure  33) 


j-  -  1  =  C 
n 


F 
(In  )2 
l  M»  ' 


7 


(37) 


If  h>*  is  substituted  for  Y,  -  Y  ,  equation  (37)  can  be  written  in 
the  following  form: 


V 


=  c 


(  —  r 


X 


(38) 


This  is  the  presentation  used  in  the  graphs  at  the  end  of  this  section. 
Figures  55  to  6t*-  present  the  generalized  backwater  equations  for  the 
geometries  considered.  Significant  values  of  the  geometries  are 

used  as  parameters. 

V      Fn  2 
Values  for  =-  and  (  r|7  )  were  calculated  in  the  same  computer 

'n 

program  as  used  for  the  head  loss  considerations  in  the  previous 

section.  These  points  were  to  be  plotted  on  a  log  log  scale,  so  the 

computer  program  was  also  written  to  determine  the  slope  of  the  best 

fitting  line  through  these  points.  A  least  mean  square  routine  for 
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log  log  plots  was  used  (see  Appendix  C).  For  the  smooth  boundaries 
a  sufficient  number  of  data  was  available  to  use  the  least  square 
method  in  finding  the  slope  of  the  lines.  In  the  case  of  rough 
boundaries,  however,  there  was  sometimes  necessary  to  fit  the  lines 
by  eye  from  the  few  data  available. 

The  data  used  for  the  calculations  are  listed  in  runs  101  -  306, 
presented  in  Tables  XVI  -  XXVII.  The  calculated  values  are  listed  in 
Tables  XXX  and  XXXI. 

There  are  some  advantages  in  doing  all  these  calculations  in 
one  computer  program.  Only  one  set  of  data  has  to  be  prepared  and 
the  computer  stores  the  values  in  the  memory  once.  The  computer  time 
used,  therefore,  becomes  a  minimum.  For  all  the  calculations  per- 
taining to  the  figures  (k6   -  65)    the  time  used  by  the  IBM  7090  was 
less  than  10  minutes. 
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FIG.    61      SUMMARY   OF    BACKWATER      RATIO      CURVES      FOR 
GEOMETRIES     la    AND   lb    .SMOOTH    BOUNDARIES 
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CHAPTER  VIII 
FINDINGS,  CONCLUSIONS  AND  RECOMMENDATIONS 

The  results  of  this  report  are  applicable  to  bridges  with  sub- 
merged inlet  and  geometries  I-a,  I-b,  V-b,  VI  and  VII  (see  definition 
of  geometries  in  Chapter  V  and  figure  2^).  No  eccentricity,  entrance 
rounding  or  wing  walls  are  counted  for. 

Findings 

1 .  Equation  of  discharge  for  orifice  type  flow  through  a 
semi-circular  submerged  bridge  constriction. 

2.  Semi-empirical  expression  for  the  occurrence  of  free 
surface  flow  and  orifice  type  flow  through  the  bridge  opening. 

3.  Graphs  presenting  the  discharge  as  a  function  of  the 
ratio  of  the  backwater  depth  to  the  arch  span. 

k.     Graphs  showing  the  ratio  of  the  backwater  to  the  normal 
depth  as  a  function  of  the  normal  depth  Froude  number  and  the  channel 
opening  ratio. 

5.  Empirical  generalized  backwater  equations  derived  for  the 
several  geometries  from  the  graphs  in  point  4, 

6.  Graphs  giving  the  value  of  the  constriction  head  loss 
coefficient  as  a  function  of  the  channel  opening  ratio  for  the 
different  test  geometries. 


131 


Conclusions 

7a.  The  backwater  at  submerged  arch  bridge  constrictions 
perpendicular  to  the  stream  decreases  as  the  length  of  the  bridge  in 
the  direction  of  flow  increases. 

7b.  For  equal  bridge  length  the  backwater  due  to  submerged 
arch  bridges  is  further  decreased  when  the  constriction  is  skewed 
with  respect  to  the  direction  of  flow,  at  least  for  skew  angles  up 
to  30°. 

8.  The  value  of  the  head  loss  coefficient  due  to  a  submerged 
arch  bridge  constriction  is  independent  of  the  normal  depth  Froude 
number  as  long  as  both  inlet  and  outlet  of  the  bridge  are  submerged. 

9»     ?or  hydraulically  short  submerged  arch  bridges  (that  is 
where  the  barrel  is  not  sufficiently  long  to  allow  the  expanding 
depth  of  flow  below  the  contraction  to  raise  and  fill  the  barrel), 
the  backwater  depth  decreases  for  an  increase  in  the  normal  depth 
Froude  number. 

10.  To  obtain  a  minimum  backwater  depth  at  a  submerged  arch 
bridge  constriction  the  barrel  should  flow  full  (full  conduit  flow). 

Recommendations 

11.  Further  experiments  as  to  the  value  of  the  head  loss  co- 
efficient through  different  types  of  bridge  constrictions  including 
eccentricity,  wing  walls,  entrance  rounding  and  roughening  of  the 
inlet  section. 

12.  Investigation  of  the  spearation  and  vortex  phenomena  at 
the  bridge  inlet. 
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APPENDIX  A 

t 

NOTATIONS 

DEFINITION 

Area 

Total  depth  water  area  at  section  1 , 

A1  =BY1 

Normal  depth  water  area,  A  =  BY 

m   n 

Total  depth  flow  area  at  section  2 
Area  of  bridge  opening 

Rectangular  channel  width 

Diameter  of  circle  of  which  the  arch 

segment  is  a  part 

b         L  Span  width  at  the  springline  of  the  arch 

C  A  coefficient 

C         LT/T         The  Chezy  roughness  coefficients  V  // HnS 

G  Coefficient  of  contraction 

c 

C,  Coefficient  of  discharge 

a 

C  Coefficient  of  velocity 

d         L  Diameter  of  a  circle 

D         L  Hydraulic  depth  as  defined  by  V.  T.  Chow 

d'         L  Distance  from  the  springline  to  the  center 

of  curvature  of  the  arch 


SYMBOL 

DIMENSIONS 

A 

L2 

A1 

I2 

An1 

I2 

An2 

L2 

A 

0 

L2 

B 

L 

b1 

L 
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SYMBOL  DIMENSIONS           DEFINITION 

E<  ^  Energy  loss  between  section  1  and  b 

e  eccentricity 

F  Denotes  a  mathematical  function 

F  Normal  depth  Froude  number  =  V  /^/gY 

f  Denotes  a  mathematical  function 

G  Denotes  a  mathematical  function 

g  Denotes  a  mathematical  function 

g  L/T         Acceleration  of  gravity 

H  L                           Total  energy  head 

h  *  L                           Y.   -  I 

1  1          n 

h  Static  head 

hf  L           Headloss  in  constriction 

K  Weir  coefficient  as  defined  in  Francis 

weir  formula 

K  Head  loss  due  to  bridge  constriction 

alone 

L  L           Length  of  the  bridge  in  the  direction  of 

flow 

L,  L           Distance  between  bridges 

L..  k  L           Distance  between  section  1  and  4 


APPENDIX  B 

DEVELOPMENT  0?  THEORETICAL  DISCHARGE  EQUATION  FOR  ORIFICE 
TYPE  FLOW  THROUGH  SEMICIRCULAR  OPENINGS 
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SYMBOL     DIMENSIONS  DEFINITION 

M  Channel  width  ratio,  b/B 

M'  Channel  opening  ratio,  A  /A 

Mi  Mild  slope  backwater  curve  in  an  open 

channel 

N  Number  of  spans 

n  L  Manning's  roughness  coefficient 

n  Subscript  which  refers  to  the  normal 

depth  for  uniform  flow 


?/L  Pressure 


Q         L-^/T        Actual  discharge  through  opening 

Q.  Theoretical  calculated  discharge  through 

opening 
Q  L3/T         Total  flow 

q         l?jl  That  portion  of  the  total  flow  which 

could  pass  through  the  bridge  without 

contraction 

R         L  Hydraulic  radius 

R  Reynolds  number  VR  /,, 

n  n  u 

r         L  Radius  of  curvature  of  the  arch 
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SYMBOL     DIMENSIONS  DEFINITION 


Slope 


S  Slope  for  normal  depth  flow  in  unconstrict- 

ed  channel 


Time 


V  L/T         Average  velocity 

V  L/T         Average  velocity  through  constriction, 


Y 

L 

Y1 

L 

Yn 

L 

Yt 

L 

a 

P 

7 

F/L3 

V 

L/T 

P 

ftV 

q/a 


o 

V  L/T         Average  velocity  when  uniform,  normal 

flow  occurs,  Q/BY 

V  Volt 

Depth  of  flow 

Depth  of  flow  at  section  of  maximum 

backwater 

Depth  of  the  normal  unconstricted  flow 

tailgate  height 

Kinetic  energy  coefficient 
Segment  factor  d'/r 
Specific  weight  of  water 
Kinematic  viscosity  of  the  fluid 
Fluid  mass  density 
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APPENDIX  3 

DEVELOPMENT  OF  THEORETICAL  DISCHARGE  EQUATION  FOR  ORIFICE 
TYPE  FLOW  THROUGH  SEMICIRCULAR  OPENINGS 

Development  of  theoretical  discharge  equation  for  orifice  type 
flow  through  Semicircular  openings.     With  reference  to  Figure  5t 
and  neglecting  the  velocity  of  approach  the  discharge  is: 


Q 


t  =  /vdA  (8) 


where 


V  =        V2g   (Y,   -  h) 


and 

h=r 


A=2  f         ^/  r2  -  h2  dh 

h=o,/ 


Hence , 


Qt  =  2       J  \/2g     (Y1    -  h)  ^/?  -  h2         dh 


=  2  V2i~       [T       (t,    -  h)»    (  r2  -  h2)^     dh  (9) 


Take  Y1  z  and  r  outside  the  integral 


Q,  =  2       v^    Y^  r     f     (  1   -  §-  )*      (  1   -  ~     )*     dh 
o  1  r 


where 
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Expanding  each  of  the  terras  in  the  integrand  into  a  binomial  Series: 
Qt  =  2       Tig"   Y^  r     /      (  1  -  i  ^  -  1/8  2L,2  -  1/16  X,3 

-  5/128  Xj4  ...)   (  1   -  i  X2  -  1/8  X22  -  1/16  x  2  -  5/128  X^..)  dh 

^   =  h/Y., 
X,  =  h2/r2 

Multiply  the  two  series  together  using  the  five  leading  terms, 
the  integrand  becomes: 

r   (  1  -  i  x1  -  1/8  x.,2  -  1/16  x13  -  5/128  x^  -  i  Xg  +  i  x1  x2 

o 

-  1/16  Xj2  ^  +  1/32  X,3  X2  +  5/256  Xj      X2  -   1/8  X22 


+  1/16  :^  x22  +  1/64  x12  x22  +  1/128  x.,2  x23  +  1/256  x^  x23 

+  5/2048  X,*  X22  -  5/128  I2k   +  5/256  X1  X24  +  5/1024  X,2  X24 
+  5/2048  X^  X2^  +  5/16,38^  X,4  l£~)     dh 


Adding  terms  and  substituting  for  X.  and  X2  in  the  integrand 


/ 


-  5/128   (  fe-  )k  +  5/256   (  fe-  )4  (  &*  )  +  5/1024  (  |U  )* 


"1 


1 


2  '    '   ■"  ,V"*T    v  Y, 
r  1 


(  £  )2  +  5/2o4<5  (  |-  )4  (^  )3  -  5/128  (  ^  )4  +  5/256 
r  1  r  r 
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C  fc-  >    (  4  >*  +  5/1024   (  £-  )2  (  ^  )^  +  5/2048   (  §-  )3 
y1  r  Y1  r2  Y1 


(  4  ^  +  5/16,384  (  *L-  )*  (  iil  )* 

r  1  r 


dh 


Integrate  and  substitute  for  limits  of  integrand. 


5r 


(  -  1/128     1/5  +  1/256     1/7  +  1/1024     1/9  ♦  1/2048     1/11 
+  1/16,384     1/13   )    (  |-  )4  +  1/2048     1/12   (  ~  )3  +  1/1024 


1/11    (  f-  )2  +  1/256     1/10      (£-)-  1/128     1/9 
*1  x1 


Simplifies  to 


-  0.00435  +  0.00195  (  7-  )  +  0.000455  (  $-  )2 
x1  x1 


+  0.00020 


(  |-  )3  -  0.00430   (  £-  )4- 


L1 


1 


Now, making  the  first  term  in  the  new  series  equal  to  unity  and 
substituting  the  solved  integrand  back  into  the  expression  for  the 
theoretical  discharge,  remembering  that     2  r  =  b 


(^  =  0.4019        V2g"    Tt*  b2    [l    -  0.2136   (  —  )  -  0.03216   (  |-  ) 


-  0.0112   (  7-  )3     -  O.OO5344   (  £-  )*...] 


1 
(10) 
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This  result,   using  the  five  leading  terms  in  the  binomial  series, 
can  be  written 

Qj.   =  C1    Y^  b2  T  (11) 

where 


and 


C1    =  0.4019  V2g     =     3.22 

T  =   1    -  0.2136    (   f-  )    "   0.03216    (   j-   )2  -   0.0112   (   j-  )3 


0.0053^*    (   7-  )4 
X1 


The  above  theoretical  discharge  equation  is  valid  for  Y..>r 
which  means  as  long  as  the  constriction  is  submerged. 
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APPENDIX  C 
DIGITAL  COMPUTERS  AND  PROGRAM  NOTES 

LGP-30  Digital  Computer 
The  LGP-30  is  a  stored  program  general  -  purpose  computer  marketed 
by  the  Royal  McBee  Corporation.  It  is  designed  for  every  day  engineer- 
ing problems  and  belongs  to  the  more  inexpensive  class  of  computers. 
The  three  basic  components  of  a  computer  memory  unit,  arithmetic 
unit,  and  control  unit  are  assembled  in  one  desk-size  housing.  The 
main  memory  of  4096  words  is  a  magnetic  drum  which  rotates  at  4000 
revolutions  per  minute,  thus  making  possible  a  maximum  access  time  of 

15  milliseconds.  The  memory  is  arranged  on  the  drum  on  64  tracks, 
each  track  containing  64  sectors.  Punched  paper  tape  is  the  permanent 
file  of  information  for  the  LGP-30.  The  primary  input-output  device 
is  a  tape  typewriter.  The  input  information  to  the  LGP-30  is  trans- 
mitted by  means  of  punched  tape  that  is  fed  into  the  computer's 
reader.  Manually  input  is  also  possible  by  writing  the  instructions 
on  a  code  typewriter,  Flexowriter,  on  the  machine. 

The  final  answer  is  printed  out  on  Flexowriter  paper.  The 
instructions  needed  to  operate  the  LGP-30  consists  of  a  series  of 

16  elementary  instructions  which  are  read  and  obeyed  by  the  computer 
one  after  another.  Instructions  are  the  orders  which  tell  the  computer 
what  we  want  it  to  do  with  the  numbers  which  we  have  given  it. 
Instructions  are  stored  inside  the  machine  before  the  operation  of 
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the  problem  is  started.  Then  these  instructions  are  considered  by 
the  machine  one  by  one  and  the  data  are  processed  accordingly.  In- 
structions have  two  major  portions.  The  first  portion  is  a  code, 
called  the  operation  code,  which  tells  the  computer  which  operation 
is  to  be  performed  next,  like  "add".  The  second  portion  tells  the 
computer  where  to  find  the  number  to  be  used  in  this  operation.  This 
portion  is  called  the  address  portion. 

Internal  calculation  is  accomplished  by  directing  to  the  computing 
section  from  memory, processing  it,  and  directing  it  back  to  memory. 
The  computing  section  consists  of  three  working  registers:  the 
counter  register,  the  instruction  register  and  the  accumulator  register. 
The  three  registers  are  located  on  three  separate  recirculating  tracks 
on  the  L3P-30  drum,  in  addition  to  the  64  tracks  of  memory. 

The  accumulator  is  the  working  register.  It  contains  the  results 
of  addition,  subtraction,  multiplication  etc.,  as  these  operations 
occur.  It  also  contains  one  of  the  operands,  prior  to  execution  of 
an  arithmetic  instruction;  the  second  operand  is  in  memory  at  the 
work  location  specified  by  the  adress  part  of  the  instruction. 

The  instruction  register  contains  the  instruction  being  executed. 
Both  the  operator  and  the  address  of  the  instruction  are  located  in 
this  register.  Words  are  transmitted  from  memory  to  this  register. 
Generally,  these  words  are  instructions  of  a  program  stored  in  memory. 
The  instruction  register  may  be  interrogated  to  determine  the  opera- 
tion to  be  performed. 

The  counter  register  contains  the  location  of  the  next  instruc- 
tion to  be  executed,  i.e.,  if  the  counter  register  reads  4132,  the 
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next  instruction  to  be  executed  is  stored  at  4132  or  track  41,  sector 
32.  This  register  contains  only  the  address  part  of  a  word.  The 
address  part  of  the  instruction  register  may  be  transferred  to  the 
counter  register.  The  address  part  of  the  counter  register  may  be 
transferred  through  a  plus  1  adder  to  memory  by  use  of  the  instruction 
register  command.  During  the  execution  of  every  operation  1  is  added 
to  the  counter  register  automatical 1 y,  so  it  will  be  ready  to  search 
for  the  next  instruction,  after  it  finishes  the  one  being  executed. 

The  program  used  in  data  calculations  here  was  written  as  instructed 
in  a  manual  published  by  the  School  of  Electrical  Engineering  at  Purdue 
University  in  April  1960.  The  manual  is  entitled,  "LGP-30  Floating 
Point  Interpretive  System  (H1-24.3)  -  Programming  Notes." 

The  machine  time  required  to  execute  and  print  out  the  answers 
for  the  data  in  the  following  program  was  about  4  hours.  The  same 
program  could  have  been  written  much  easier  on  the  I EM  1620  or 
Fortran  7090  and  the  operation  time  would  only  have  been  a  percentage 
of  the  time  the  LGP-30  used.  The  reason  for  using  the  LGP-30  was  to 
compare  the  efficiency  of  the  LGP-30  to  the  much  more  expensive  IBM 
machines  mentioned.  Hence  for  all  other  programs  the  IBM  machines 
were  used. 
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Program  for  Calculation  of  Normal  Depth  when  Discharge.  Channel  Width 
and  Slope,  and  Mannings  Roughness  Coefficient  are  Known 

The  mean  velocity  of  a  free  stream  in  uniform  steady  state  condi- 
tion can  be  computed  by  the  Manning  formula: 


=  U^86r2/3s1/2  (3?a) 


V  = 


where:  V  =  mean  velocity  in  fps  (V  =  Q/A) 

n  =  roughness  coefficient 

R  =  hydraulic  radius  in  f t  (R  =  A/P) 

S  =  slope  of  energy  line 

Q  =  discharge  in  cfs 

A  =  cross  sectional  area  in  sq.  ft. 

P  =  wetted  perineter  in  ft. 
when  applying  this  formula  to  a  rectangular  channel: 

A  =  By    and    P  =  E  +  2y 

where  d  =  width  of  channel  in  ft. 
y  =  depth  of  flow  in  ft. 

when  uniform  steady  state  flow  occurs  y  =  y 

Substituting  into  equation  (37a)  and  rearranging  gives 

^n  =  -573  -^772  (  ^T^  ^  (37b) 

n   B5/3  1.486S1'2      yn 


The  formula  can  not  be  solved  for  y  explicitly  so  successive 
approximation  was  used.  An  approximate  depth  y„  was  chosen  and 
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equation  (37b)  which  is  of  from  Y  =  f(y)  was  calculated  and  y,  =  f(y) 
was  obtained.  Then  by  doing  the  computations  with  the  new  yf  a  y? 
value  was  obtained  etc.  until  the  change  of  y  in  percentage  is  equal 
to  or  less  than  a  given  €   . 

^  -  ym-1 


<  € 


Vi 


As  a  first  approximation  for  the  depth,  y  =  B  was  used. 
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Location 

4000 
4001 
4002 
4003 
4004 
4005 
4006 

4007 
4008 
4009 
4010 
4011 
4012 

4013 
4014 
4015 
4016 
4017 
4018 

4019 


Prop-ram 

Address 

Remarks 

; 0004000 

r6300" 

Special  Instruction 

U04001 

n                    n 

i5000' 

b5000' 

S 

rOOOO1 

Vs~ 

h5018' 

Store  vS 

b500 2' 

B 

e5010' 

35/3 

r5000' 

Q/B^/3 

m5004' 

Q/3^3 

d5012" 

Qn/ 1.486  B5^3 

d5018« 

Qn/1.486  35/3  S1/2 

h5020' 

Store  Q  / 1.486  B^3  S1/2 
n' 

b5002' 

3 

h5022' 

Store  B  as  yn 

O5022' 

yo 

h5024' 

Store  yQ 

m50l6' 

2^o 

a5002« 

B  +  2yQ 

d5022' 

B  +  2yQ 

4020 


e50l4« 


B  +  2y^ 


i/3 
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Location 


4021 


4022 
4023 

4024 
4025 
4026 

4027 

4028 
4029 
4030 
4031 
4032 
4033 
4034 
4035 
4036 
4037 
4038 


Address 
IT150201 


h5022' 
s5024« 

d5024' 
tOOOO' 
a5O08' 

t4015* 

P50001 
dOOOO1 
p5002» 
dOOOO1 
p5004' 
dOOOO' 
p5006' 
dOOOO' 
P5022' 
mOOOO ' 
u4002» 

.000  4002' 


y^  ■ 


Remarks 
Qn 


1   "  1.486  B^  SV2 


11^2)2/3 


Store  y1 


y1    "  y0 

yi  "yo 


■»( 

3 

'l  "  yo 

- 

yo 

y1  "  yo 

t 

If 

yo 
h  -yo 

If 

yo 

Print  S         * 

Tab 

Print  3 

Tab 

Print  n 

Tab 

Print  Q 

Tab 

Print  y 

Car.  Ret. 

Start  with  new  set 
of  data 

Constant  data 
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; 0004002'  Instructions  for  data  tape  to  move  up  and 

i5006'  give  new  Q  value  after  each  comouted  Y  value 

i-     n 

; 0004023' 

u4034'   Go  to 
; 0004034' 

z5006'  Print  Q  from  memory 

dOOOO'  Tab 

z5022*  Print  computed  Yn 
; 000401 V 

D5028'  B/10  (instead  of  5002) 
;  0005040' 

r6300'  Special  instruction 

u0400' 

15028'  Input  to  memory  (B/10) 

u4002'  Go  to 
; 0005040 ' 
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Data  Location 


Location 


Content 


Constant 


Variable 
Constant 


5000 
5002 
5004 
5006 
5008 
5010 
5012 
5014 
5016 
501 8 
5028 


Slope  S  =  0.001 

Channel  width  B  =  4.955  ft. 

Mannings  n  =  0.0110 

Flow  Rate  Q 

Tolerance  e  =  0.001 

5/3 

1.486 

2/3 
2 

</s  =  0.25 

B/10 
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The  IBM  1620  Computer  Using  Fortran  Programming 
The  word  FORTRAN  is  derived  from  FORmula  TRANslation.  It  is 
provided  by  I EM  for  the  IBM  1620  Data  Processing  System  and  it  allows 
engineers  and  scientists  to  write  their  problems  and  formulas  in  a 
language  closely  akin  to  their  normal  mathematical  notation.  The 
IBM  1620  then  translates  this  language  into  a  series  of  instructions 
which  are  used  to  solve  the  problem.  Because  the  FORTRAN  system 
assumes  the  responsibility  of  translating  the  problem  for  computer 
solution,  the  user  of  the  IBM  1620  need  not  be  exessively  concerned 
with  machine  details. 

Problems  can  be  organized  and  specified  for  machine  solution  more 
quickly  in  FORTRAN  language  than  in  machine  instructions  such  as  the 
LBP-30  in  using.  The  FORTRAN  language  is  written  as  a  compromise 
between  the  language  of  the  computer  and  the  language  of  the  engineer 
and  scientist.  To  satisfy  the  engineer  and  scientist,  as  many  of  the 
detailed  computer  control  operations  as  possible  are  eliminated  from 
the  job  of  writing  programs,  and  a  problem  statement  format  close  to 
that  of  the  mathematical  equation  is  used.  This  condensed  symbolic 
and  mathematical  form  of  the  problem  reduces  the  possibility  of 
clerical  errors,  highlights  and  program  logic,  and  makes  it  possible 
to  quickly  modify  or  change  the  problem. 

A  problem  written  in  the  FORTRAN  language  if  made  up  of  state- 
ments which  can  be  identified  by  an  arbitrary  number  if  so  desired. 
Three  main  types  of  statements  are  (1 )  arithmetic,  (2)  input-output, 
and  (3)  control.  Arithmetic  statements  provide  for  computation  of 
the  problem.  An  example  is  C  =  A  +  B.  This  statement  would  cause 
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IBM  1620  to  generate  a  series  of  machine  instructions  which  would 
set  the  variable  C  equal  to  the  sum  of  the  variables  A  and  B.  Other 
types  of  arithmetic  operations  are  multiply,  represented  by  *,  divi- 
sion, represented  by  /,  and  exponentiation  represented  by  **  (e.g. 
A**2  is  the  same  as  A  ).  The  use  of  these  operations  would  be  illustra- 
ted by  the  FORTRAN  statement: 
X  =  A*B  +  C/Y  -  Y**3 
which  is  equivalent  to  the  mathematical  expression 1 

X  -  AB  +  C/Y-Y3 
Input  -  output  statements  provide  for  transfers  of  data  into  and 
results  out  of  the  computer.  A  type  of  output  statement  is  PRINT 

I ,  X,  Y  which  would  cause  the  values  of  X  and  Y  to  be  printed.  Control 
statements  regulate  the  sequence  in  which  calculations  are  performed. 
Two  of  the  control  statements  are  STOP  and  IF.  The  STOP  signifies 

the  end  of  the  problem.  The  IF  statement  provides  for  program 
decision  depending  upon  whether  or  not  the  specified  variable  is  less 
than,  equal  to,  or  greater  than  zero.  The  statement  IF  (X  -  Y)  10, 

II,  12  will  cause  instructions  to  be  compiled  that  will  test  the  value 
of  X-  I,  If  the  value  is  less  than  zero  (X<Y),  statement  10  will 

be  performed  next.  If  X  -  Y  equals  zero  (X  =  Y),  then  statement  11 
will  be  performed  and  if  X  -  Y  is  greater  than  zero  (X>Y),  then 
statement  12  will  be  performed. 

Cost  Comparisons  -  IBM  1620,  LGP-30,  IBM  7090 
The  IBM  1620  is  a  solid  state  digital  computer,  which  means 
that  the  computer  is  built  with  primarily  solid  state  devices  such 
as  transistors,  distinguished  from  non-solid  state  devices  such  as 
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vacuum  tubes  which  is  the  main  structure  of  the  LGP-30.  Solid  state 
devices  are  generally  more  reliable  than  non-solid  state  devices. 
The  I EM  also  has  an  automatic  built-in  checking  program  that  detects 
and  permit  corrections  of  the  more  common  programming  errorso  LGP-30 
has  no  automatic  checking. 

The  cost  of  the  LGP-30  is  about  $50,000  while  the  IEM  sells  for 
from  $75,000  and  up  to  $200,000. 

Another  advantage  of  the  IBM  1620  is  that  a  problem  specified 
in  the  FORTRAN  language  for  this  machine  can  be  run  with  modifications 
on  other  IBM  data  processing  systems  using  FORTRAN  systems.  Other 
programming  systems  such  as  GOTRAN,  Symbolic  Programming  System  and 
package  subroutines  can  also  be  used  by  the  IBM  computers.  Other 
IBM  computers  such  as  the  IBM  7090,  also  used  for  data  computation 
in  this  study,  costs  about  $3,000,000.  IBM  7090  is  also  a  solid  state 
digital  computer,  but  its  price,  speed  of  computation  and  capacity 
of  memory  makes  it  more  suitable  for  a  larger  concern  or  Research 
Institute. 

The  following  programs  are  written  in  the  FORTRAN  language  and 
processed  on  the  IBM  1620  or  7090  computers. 

Program  for  Prandtl  Tube  Calibration 

Velocity  =  f  (Voltage) 
Assuming  the  relation: 

Ah  =  BV 

A  h  is  in  cm. 

Let:       Ah  =  Y 
V  =  X 
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Then  Y  =  B  X 

Now  compute  B  by  least  square  method: 

2  £  (Y,  -  BX,)2 

£  (Y,   -  BX,)2  =  min  -     S_£ i !L_    =     0 

dB 


Y,  2   (Yi  -  BXi)    (  -  Xi   )  =  0 

VX.    Y.    -  B  VX  2  =  0 

^ii  *-•    i 


L     X,   Y, 

B  "      n       T  2 


The  correlation  coefficient 

n 

R 


n        1     n        n 

I  X.  Y.  -  -    T    X.     V  Y4 
Viin    x_i    x_i 


xy 


[  S\2-  £<    ?\>2]    [  i \2-£<  IV2] 


Since: 


Vy  =  C  V2gAh 


where : 


Vv  is  in  fps. 


Ah  is  in  feet 
i.e.   for  Prandtl  tube: 


Vv  =  V2gAhft 
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converge   A  h  from  ft  to  cm : 


V/^3  =V//2g    (  0.032808)  A  hran 
Vv  =     /2  x  32.174  x  0.032808  3V 

v   =  c/v" 

v 


where 


C  =*Jz  x  32.174  x  0.032808  B 
The  program  is  written  to  compute  E,   R%  C  for  the  measured  data. 
Program     IBM  1620 

2  DATA  FITTING  FOR  VELOCITY  =  F(VOLTAGS)  BY  Iffl  1620. 

DIMENSION  X(300),Y(300) 
1 7  READ  1 , N , J 

1  FORMAT   (15,15) 
TYPE   2 

2  FORMAT (//////24H  H(CM)  =  B  TIMES  VOLTAGE) 
TYPE  3 

3  FORMAT  (48H  VELOCITY  (FPS)  =  C  TIMES  SQUARE  ROOT  OF  VOLTAGE//) 

TYPE  4 

» 

4  FORMAT  (3X,1HJ,10X,1HB,11X,3HRXY,11X.1HC/) 
CONST  =  2.0*32.17^*0.03280033 

K  =  1 
L  =  J 
M  =  1 
P  =  J 
14  SXY  =  0.0 
SXSQ  =0.0 
SX  =  0.0 
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SY  =0.0 
SYSQ  =  0.0 

10  DO  7  I  =  K,J 

SXY  =  SXY  +  X(I)*Y(I) 

SXSQ  =  SXSQ  +  X(I)**2 

SX  =  SX  +  X(I) 

SY  =  SY  +  Y(I) 
7  SYSQ  =  SYSQ  +  Y(I)**2 

B  =  SXY  /  SXSQ 

RXY  =  (  SXY-SX*SY/P)/SQRT((SYSQ-SY**2/P)*(SXSQ-SX**2/P)) 

C  =  SORT (CONST  *  3) 

TYPE  8,J,B,RXY,C 
3  FORMAT  (15.F13.5.F13.5.F13.5) 

IF  (J  -  N)  9.11,15 
9  K  =  J  +  1 

J  =  J  +  1 

P  =  J 

GO  TO  10 

11  00  TO  (12,17),M 

12  TYPE  13 

13  FORMAT  (/) 
K  =  L  +  1 
P  =  J  -  L 
M  =  2 

GO  TO  14 
15  TYPE  16 
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16  FORMAT   (6H  2RR0R) 
GO  TO  17 
END 

Program  for  calculation  of  curves  in  test  selection  charts,  (see  fig.  26  &  27) 
Development  of  tailgate  hight  eqn.  from  Francis  weir  formula: 

(37) 


Q  =  KB 


(H+||)3/2.(  ^/Z 


where  Q  =  discharge  cfs. 

B  =  effective  length  of  crest,  ft. 

H  =  measured  head  on  crest,  ft,  upstream  from  weir 

beyond  beginning  of  surface  curve. 
V  =  velocity  of  approach 
K  =  weir  coefficient 


Set 


Y,  =  tailgate  hight 


Y  =  normal  depth  of  flow 


Then  eqn  (37)  can  be  written  in  the  following  way: 

£-=   (  Y     -  Y     +£bv/2  .    (  2&)3/2 
KB        k     n       h       2g   ;  {   2g   ; 

rewrite   (38) 


(36) 


S.+  (  2^)3/2  ] 

L  KB        l  2g   ;         J 


2/3 


V2 


■     Yn"Yt  +  2ia 


Now  solve  for  Y, 


Yt=Yn  +  Qii     " 


h    +   <4  '^  ] 


2A 
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but 


and 


V  = 


A   3Y 
n 


y  =  a 


2g  B' 


■3 +  Y  - 

2  _  2      n 


a.    +  (  Q^  )3/2 

KB      '  2g  B2  Y  2 

n 


2/: 


'(39) 


when  substitution  the  constant  values 
Ct  =  assumed  to  be  unity 
K  =  ^. 205 
3  =  4,955  ft. 

The  tailgate  equation  is  found  to  be 


\  = 


1581   Y 


-x       +  Y 
2  n 


20^     +      < 


1581    Y 


,3/S 


n2/3 


(40) 


The  Y.    values  are  calculated  in  the  urogram  and  dotted  versus  Y 
t  i      t>  -  n 

with  Q  as  parameter  in  the  nomographic  charts  in  fig.  26  and  2?. 


Program     IBM  162 


20  DO  10  1=10,45,5 
PN=0.0 
PI=I 
G=.1*PI 

PRINT  4 
4  FORMAT  (30H 
8  READ  6,Y 
6  FORMAT  (F4.0) 

PN=PN+1.0 


YN 


YT 
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N=PN 
D=1. 
A=(D*(Q**2/((1581.0)*(Y**2)))+Y) 

&=((D*(Q**2/(1531.0*(Y**2))))**1.5) 
C=((Q/ 20.836)+E)**(2./3.) 
YT=A-C 

PUNCH  2.Q.Y.YT 
2  FORMAT   (F10.5,2X,F10.5,2X,F11.6) 
IF   (N-42)   8,10,10 
10  CONTINUE 
PAUSE 
GO  TO  8 
END 

Program  to  find  Theoretical  Discharge  Q, ,   and  Coefficient  of  Discharge  C,« 

From  equation  (11),    Chapter  III 

0^  =  C1   Y^  b2T  (11) 

where 

C1    =  3.22  Cd 


and 


T  =  1    -  0.1068   (  ~  )  -  0.00804  (  ^-  )2  -  0.0014  (  £■  )3 
X1  *1  Z1 

-  0.000335  (  7-  )4 


Mote:  b  =  2r 
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Equation  (11)  is  programmed  and  with  b,  Y, ,  Y  ,  Q  and  F  as 
given  data,  the  print-out  will  be: 

V  Cd'  <W<VYn^  Fn»  W  ^^ 

Program     1 5-1  1620 

PRINT  80 
80  F0RMATO9H  Q?     CD     A0/AN*YI/YN     FN     YI/YN     (FN*AN/A0)**2 
20  READ  10,QA,F,YN,Z,P 
10  FORMAT  (F2.1,2X,F2.2,2XtF4.2,2XtF3.1,2XtF4.2) 

B=4.955 

YN=YN*(.0328) 

H=((o0328)x(P)) 

H=((.0328)*(P)) 

Y=((.0833)*(Z)) 

C=((. 000335)* ((Y/H)**4)) 

A=(1.-(.1068*(Y/H))-(.008o4*((Y/H)**2))-(.0014*(Y/H)**3))-C 

D=((3.22)*((H)**(.5))*(Y**2)) 

QT=A*D 

cd=qa/qt 

AN=YN*B 

A0=(3.14/8.)*(Y**2) 

AI=H*B 

aa=(a0/an)*((h/yn)**.5) 
ab=h/yn 

AC=((F*AN)/A0)**2 
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PRINT  90,QT,CD,AA,F,AB,AC 
90  FORMAT   (F4.1,F7.3,F8.4,F1^.2,F7.3,F9.3) 

IF  (Z-18.)   20,100,100 
100  PRINT  110 
110  F0R!-IAT(49H ) 

GO  TO  20 

3ND 

Program  for  calculation  of  the  plotting  values  for  the  log  log  graph 

F 
with  axis  (  =  -  1  )  versus  (  rr  )  ,  and  C . 
n 

Program  IBM  1620 

PRINT  80 
80  F0RMiAT(^9H  QA     CD     AC/AN*(YI/YN)£   FN     Yl/YN        (FN*AN/AC) 
20  R2AD  10,QA,F,Y!I,Z,P 
10   FORMAT   (F2.1,2X,F2.2,2X,F^.2,2XtF3.1,2X,F4.2) 

B=4.955 

YN=YN*(.0328) 

H=((.0328)*(P)) 

Y=((.0833)*(Z)) 

C=((.OO0335)*((Y/H)**if)) 

A=(1.-(.1o68*(Y/H))-(.0080^*((Y/H)**2))-(.001^*(Y/H)**3))-C 

D=((3.22)*((H)**(.5))*(Y**2)) 

QT=A*D 

CD=QA/QT 

AN=YN*B 

AC   (3.1V8.)*(Y**2) 
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AI=H*B 

aa=(ac/aii)*((h/yn)**.5) 

AB=H/YN 
AC=((F*AN)/A0)**2 

PRINT  90,QA,CD,AA,F,AB,AC 
90  FORMAT  (F4.1,F7.3,F8.4,F1^.2,F7.3,F9.3)' 

IF  (Z-18.)    20,100,100 
100  PRINT  110 
110  FORMAT   (4QH . — . 

GO  TO  20 

EI© 


Program  for  finding  normal  depth  Froude  number,  F  ,  when  given  discharge, 


Q,  and  normal  depth  Y  . 

Program  equation:  Fn  =  g-|-  /(g  y  V*  =  0,0^|Q 

n/  Yn^ 

Program     I5M  1620 

PRINT  20 
20  FORMAT   (49H  QA  YN  FN  ) 

40  READ  10.QA.YN 
10  FORMAT   (F8.7,2X,F8.8) 

FN=(  ( .0355*QA.)  /  (YN**1 . 5) ) 

PRINT  30,QA,YN,FN 
30  FORMAT  (F6.3,F9.5,F8.4) 

GO  TO  40 

END 
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V 
Program  for  calculation  of  plotting  values  K  versus  M* ,  —  versus 

n 

Fn  2 
(  w*  )  and  the  slope  of  the  best  fitting  line  through  the  points  on 

arithmetic  and  log  log  scale  respectively. 


The  programmed  equations  are 

A 

(35) 


V 

k  =  ~ —  -  a 

V/2g 


(^)2-(^)2 
Ay      A1 


and 


K  =  C1  +  C2  M» 


where  C.  and  C?  are  the  constants  computed  to  determine  the  slope  of 
the  line  (see  figures  k6  -   55). 

n 
and 

b.  =  C  (  f  )2X  (38) 

n 

where  C  and  y  are  the  constants  computed  to  determine  the  slope  of 

the  line  (see  figures  56-65) 

Program  IW:   70?0 
1  AKS=0.00 

1  AKS=0.00 
AMIS=0.00 
AM5S=0.00 
AKAKS=0.00 
ZELSUM=0.00 
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ZILS=0.00 

ZELSS=0.00 

ZELZS=0.00 

M=1 

READINPUTTAPES.90, ALPHA, PHI, G, B 

90  FCRMAT(4F10.5) 

100  RSADINPUTTAP3S.91 ,I,Q,YNI,YI,A0,PAR 

91  F0RMAT(15,F5.2,4F10.4) 

yn=yni/phi 

AHT=A0/(YN1*B) 
ZTA=SQRIF(G*YNI) 
FN=Q/(B*YN1*ETA) 
RH0=B*YN1 

v=q/ao 
z=2.0*g*phi*(yi-yh)/v**2 

pi=b*phi*yi 

0MSGA=A0/PI 

AK=Z- ALPHA* (AMI**2-0MEGA**2) 

AMS=AMI**2 

AKAM=AK*AMI 

ZE=(YI-YN)/TN 

ZE-(FN/AMI)**2 

ZIL=L0GF(ZI) 

ZEL=L0GF(ZE) 

ZSLS=ZEL**2 

ZELZ=ZEL*ZIL 
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WRITEOUTPUTTAPE6,92,ItAK,AMIfAMSfAKAMtFNtPHY 
VJRITE0UTPUTTAPZ6 ,  92, 1 ,  a ,  Z3 ,  2IL,  ZEL,  ZELS ,  ZELZ 

92  FOPJIAT(15.6F15.5//) 
AKS=AKS+AK 
AKCS=AMIS+AMI 
AMSS=AMSS+AMS 
AKAMS=AKAMS+AKAM 
ZILS=ZILS+ZIL 
ZELSS=ZSISS+ZELS 
ZELSUM=ZELSUM+ZEL 
ZSLZS=ZSLZS+ZELZ 
M=M+1 
IF(M-17)20,20,21 

20  GOT0100 

21  »JRITEOUTPUTTAPE6,93,M,AKS,AMIS,AMSSfAKAHS 
VIRIT30UTFUTTAPE6, 93.M,  ZIIS,  ZEISUM,  ZEISS,  ZSLZS 

93  FORMAT (1  5,^15.5//) 
CONSB=((AMIS*AKS)-(1?.0*AKAI-!S))/(AMIS**2-17.0*AMSS) 
C0NSA=(AKS-AMIS*C0NSB)/1 7.0 

CON&=(  ( ZELSUM*ZILS)=(  1 7.0*ZELZS) )  /  (ZELSUM**2-1  7.0*ZELSS) 
C0NA= ( ZIIS- ZS ISUM*  CONB ) / 1 7 . 0 
1WRITE0UTPUTTAPE6,9ii-,C0NSA,C0NS3,C0NA,C0Iffi 

94  FORMAT(4F15.5) 
G0T01 

END(1, 0,0, 0,0, 0,1, 0,0, 1,0, 0,0,0,0) 


APPENDIX  D 
TABLES 
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TABLE    XXX      CALCULATED     VALUES-   RUN     NO. 

Run  A,  Qt  Cc  Cd  C,. 


2-140 


.635  8.317 

.632  7.291 

.628  6.737 

.617  6.585 

.624  6.229 

.614  5.897 

.629  5.493 

.654  5.198 

.641  4.917 

.594  8.265 

.004  7.315 

.614  6.687 

.622  6.163 

.623  5. 450 

.515  7.232 

.48e  6.475 

.490  6.371 

.491  6.148 

.500  5.563 

.505  5.291 

.506  5.000 

.510  4.950 

.522  4.668 

.524  4.488 

.525  4.133 

.511  4.626 

.508  4.92C 

.501  5.252 

.496  5.701 

.488  0.283 

.411  6.240 

.410  5.939 

.417  5.729 

.425  5.308 

.413  5.145 

.421  4.96e 

.422  4.716 

.427  4.445 

.429  4.221 

.435  3.899 

.434  3.653 

.433  3.511 

.432  3.188 

.412  6.256 

.414  4.842 

.415  4.374 

.430  3.741 

.434  3.273 

.261  4.256 

.262  4.097 

.262  3.368 

.263  3.700 

.263  3.560 

.267  3.388 

.263  3.286 

.263  3.C88 

.263  2.964 

.263  2.792 

.264  2.706 

.273  2.b54 

.276  2.502 

.252  2.282 

.275  2.204 

.277  2.093 

.277  1.929 


2.132 
1.886 
1.955 
2.232 
2.210 
2.230 
1.960 
2.201 
3.004 
2.146 
2.020 
2.095 
1.862 
1.907 
2.188 
2.283 
2.190 
2.066 
2.807 
2.045 
2.171 
2.239 
2.040 


2.121 
2.178 
2.097 
2.040 
2.141 
2.065 
1.886 
2.211 
2.204 
1.878 
1.617 
1.555 
1.347 
2.606 
1.865 


.141 
.171 
.203 


.203 
.120 
.188 
.227 
.271 


.156 
.228 
.277 
.329 


.370 
1.062 
.758 
.515 
.844 
.594 
1.981 


1.263 
.889 
.607 
1.479 
1.063 
.727 
1.100 
1.494 
1.565 
.953 
.821 


.570 

.842 

.901 

.435 

.567 

.554 

.812 

1.120 

.393 

.583 

.552 

.798 

1.257 

.37? 

.580 

.544 

.785 

1.293 

.364 

.579 

.542 

.775 

1.400 

.352 

.578 

.540 

.771 

1.500 

.338 

.376 

.526 

.752 

1.640 

.318 

.573 

.520 

.744 

1.759 

.305 

.570 

.515 

.736 

i.e7o 

.293 

.569 

.57C 

.842 

.919 

.433 

.387 

.550 

.806 

1.110 

.392 

.583 

.556 

.800 

1.267 

.374 

.580 

.540 

.771 

1.418 

.348 

.577 

.529 

.756 

1.657 

.317 

.572 

.256 

.385 

.800 

.401 

.584 

.370 

.518 

.950 

.372 

.580 

.394 

.54b 

.980 

.367 

.580 

.440 

.613 

1.025 

.357 

.573 

.623 

.902 

1.208 

.332 

.575 

.747 

1.040 

1.297 

.321 

.573 

.917 

1.274 

1.400 

.307 

.571 

.960 

1.371 

1.420 

.306 

.571 

1.240 

1.671 

1.530 

.29! 

.568 

1.386 

1.915 

1.610 

.282 

.567 

1.949 

2.780 

1.780 

.267 

.560 

1.257 

1.793 

1.550 

,.293 

.569 

1.060 

1.435 

1.430 

.309 

.571 

.784 

1.162 

1.310 

.323 

.573 

.578 

0.843 

1.165 

.341 

.576 

.408 

0.605 

1.010 

.363 

.579 

.591 

.se3 

.690 

.371 

.580 

.589 

.879 

.750 

.358 

.578 

.595 

.8S9 

.797 

.353 

.578 

.584 

.872 

.900 

.332 

.575 

.555 
.549 
.540 
.591 
.585 
.575 


.578 
.580 
.575 
.575 
.574 
.560 


.974 
.688 
1.154 
1.277 
1.401 
1.011 
1.011 
1.011 
1.011 
.415 
.415 


.916 
.304 
.794 
.784 

.768 


.896 
.693 

.834 
.375 


ifev 


.441 
.253 
.199 

.160 


.945 
•1.000 
1.082 
1.182 
1.273 
1.420 
1.549 
1.628 

1.-:.  5 
.086 

1.040 


.650 
.695 
.755 
.792 


1.080 
1.112 
1.215 
1.380 
1.448 
1.548 
1.710 


.612 
1.160 
1.020 
.764 
1.940 
.966 
.720 
.600 
1.280 
.992 


1.960 
.982 
.769 


.820 
1.480 
.830 
1.490 
1.560 
1.455 
1.430 
1.348 
.988 
.890 


.251 
.236 


2.290 
1.7o2 
1.480 
3.330 

1.820 
1.479 
1.227 
2.308 
2.050 
1.527 

3.ee2 
1.934 
1.443 
1.201 
2.562 
1.989 
1.771 
2.300 
1.825 
3.920 
1.965 
1.542 
1.222 
3.860 
1.925 
1.471 
1.148 
2.301 
1.755 
1.367 
2.820 
2.135 
1.639 


.573 
.571 
.508 
.5o6 
.563 
.5ol 
.553 


.559 
.558 
.556 


.548 
.544 
.554 
.540 


E77 

1.395 
1.272 
1.166 
1.521 
1.219 
1.109 
1.030 
1.400 
1.274 
1.168 
1.321 
1.067 
.969 
.888 
1.242 
1.131 
1.037 
1.250 
1.14o 
1.321 
1.067 
.9o9 


1.138 
1.035 
.949 


5.404  .462 

5.789  .431 

6.417  .339 

3.874  .516 

4.047  .430 

5.124  .390 

5.522  .362 

5.447  .458 

6.397  .390 

6.517  .383 

4.226  .591 

4.854  .515 

5.038  .496 


.461 


.478 


5.428 
5.795 
6.272 
6.649 

6.389  .508 

7.375  .474 

4.246  .588 

4.906  .509 

5.282  .473 

5.510  .453 

4.214  .593 

4.840  .516 

5.106  .489 

5.225  .478 

5.434  .551 

5.771  .519 

6.010  .499 

6.136  .570 

6.565  .533 

6.852  .510 


TABLE  XXXI  -  CALCULATED  VALUES  -  RUN  NO.   HI  -  306 
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LU 

1.527 

.154 

.390 

.415 

.202 

.843 

142 

1.466 

.154 

.306 

.415 

.202 

.835 

143 

1.431 

.223 

.144 

.200 

.128 

.610 

144 

1.333 

.223 

.134 

.200 

.128 

.600 

145 

1.227 

.223 

.123 

.200 

.128 

.60C 

146 

.893 

.223 

.090 

.200 

.128 

.588 

147 

2.430 

.134 

.tee 

.5o3 

-253 

1.100 

148 

1.577 

.134 

.447 

.563 

.253 

1.065 

149 

1.503 

.134 

.427 

.563 

.253 

1.050 

150 

1.403 

.134 

.398 

.563 

.253 

1.030 

151 

1.891 

.113 

1.719 

1.804 

.441 

1.900 

152 

1.788 

.113 

1.025 

1.804 

.441 

1.839 

153 

1.684 

.113 

1.531 

1.804 

.441 

1.771 

154 

1.602 

.113 

1.457 

1.8C4 

.441 

1.718 

155 

1.941 

.17° 

.728 

.740 

.253 

.908 

156 

1.672 

.17o 

.027 

.740 

.253 

.910 

157 

1.523 

.170 

.572 

.740 

.253 

.380 

153 

1.283 

.176 

.483 

.740 

.253 

.823 

159 

1.440 

.254 

.262 

.350 

.160 

.539 

loo 

1.000 

.254 

.182 

.350 

.160 

.553 

161 

.922 

.254 

.103 

.350 

.100 

.545 

162 

.746 

.254 

.136 

.350 

.160 

.531 

163 

1.807 

.15, 

.852 

.933 

.303 

1.187 

164 

1.488 

.154 

.703 

.933 

.303 

1.091 

165 

1.474 

.154 

.696 

.933 

.3C3 

1.086 

166 

1.109 

.154 

.525 

.933 

.303 

.930 

167 

1.487 

.223 

.340 

.449 

.192 

.711 

168 

1.375 

.223 

.314 

.449 

.192 

.702 

169 

1.104 

.223 

.206 

.449 

.192 

.680 

170 

.959 

.223 

.220 

.449 

.192 

.650 

171 

1.88o 

.140 

1.100 

1.154 

.354 

1.480 

172 

1.706 

.140 

.995 

1.154 

.354 

1.400 

173 

1.500 

.140 

.875 

1.154 

.354 

1.320 

174 

1.442 

.14C 

.e42 

1.154 

.354 

1.300 

175 

1.104 

,140 

.640 

1.154 

.354 

1.160 

176 

1.499 

.202 

.424 

.550 

.224 

.835 

177 

1.455 

.202 

.411 

.55o 

.224 

.e2e 

178 

1.256 

.202 

.355 

.556 

.224 

.795 

179 

1.077 

.202 

.305 

.556 

.224 

.765 

180 

2.211 

.356 

1.494 

1.277 

.256 

.830 

181 

1.247 

.184 

.420 

.6o2 

.250 

.912 

182 

.919 

.184 

.310 

.002 

.256 

.840 

183 

1.993 

.135 

2.173 

2.10C 

.441 

1.856 

184 

2.197 

.211 

.991 

.880 

.2  53 

.930 

185 

2. 009 

.211 

.907 

.886 

.253 

.892 

186 

1.681 

.304 

.370 

.427 

.loO 

.573 

187 

1.375 

.223 

.314 

.4  50 

.160 

.557 

188 

1.555 

.100 

.790 

l.-ll 

.100 

.550 

189 

1.728 

.304 

.380 

.427 

.160 

.538 

190 

1.899 

.137 

1.092 

1.131 

.303 

1.104 

191 

2.204 

.170 

1.565 

1.401 

.354 

1.490 

192 

1.805 

.170 

1.284 

1.401 

.354 

1.324 

193 

2.202 

.100 

1.117 

1.011 

.354 

1.690 

194 

2.618 

.154 

.546 

.415 

.204 

.990 

195 

2.749 

.187 

.392 

.283 

.159 

.807 

196 

2.500 

.223 

.252 

.200 

.126 

.669 

197 

1.500 

.112 

1.352 

1.789 

.441 

I.60O 

198 

1.949 

.175 

.724 

.735 

.253 

.974 

199 

2.082 

.212 

.529 

.501 

.199 

.784 

200 

1.650 

.252 

.297 

.353 

.160 

.610 

201 

1.599 

.158 

1.312 

1.617 

.353 

1.155 

202 

1.370 

.135 

1.499 

2.159 

.441 

1.463 

203 

1.515 

.211 

.686 

.380 

.253 

.305 

204 

1.603 

.256 

.497 

.004 

.199 

.635 

205 

1.443 

.140 

.840 

1.154 

.354 

1.300 

206 

1.640 

.100 

.836 

1.011 

.353 

1.470 

207 

1.398 

.100 

.710 

1.011 

.353 

1.370 

208 

2.589 

.154 

.540 

.415 

.202 

.980 

209 

1.471 

.154 

.307 

.415 

.202 

.837 

210 

1.199 

.223 

.121 

.200 

.128 

.597 

211 

1.398 

.158 

1.150 

1.017 

.441 

1.465 

212 

1.199 

.223 

.121 

.200 

.441 

1.330 

213 

.960 

.112 

.868 

1.789 

.441 

1.320 

214 

1.939 

.175 

.720 

.734 

.253 

.970 

215 

1.336 

.175 

.498 

.734 

.253 

.845 

216 

1.119 

.175 

.417 

.734 

.253 

.801 

218 

1.419 

.211 

.643 

.386 

.354 

1.250 

219 

.941 

.140 

.551 

1.154 

.354 

1.095 

220 

.723 

.140 

.424 

1.154 

.354 

1.003 

221 

1.674 

.252 

.3d 

.353 

.160 

.013 

R-ir. 
No. 

K 

K 

hl'" 

,fr     2 

«'  b- 

7-     Jl 

b 

222 

1.530 

.252 

.285 

.353 

.253 

.605 

223 

1.085 

.252 

.196 

.353 

.160 

.562 

224 

1.398 

.158 

1.149 

1.017 

.353 

1.075 

225 

1.240 

.15e 

1.021 

1.017 

.353 

1.010 

226 

1.C79 

.158 

.890 

1.017 

.353 

.943 

227 

1.073 

.135 

1.177 

2.159 

.441 

1.270 

228 

1.036 

.135 

1.137 

2.159 

.441 

1.250 

229 

1.419 

.211 

.043 

.880 

.253 

.770 

230 

1.011 

.211 

.460 

.880 

.253 

.084 

231 

1.059 

.304 

.234 

.427 

.160 

.480 

232 

.654 

.100 

.333 

1.011 

.j53 

1.067 

233 

.251 

.223 

.025 

.200 

.128 

.504 

2  34 

.694 

.100 

.353 

1.011 

.353 

1.082 

235 

.lo9 

.223 

.017 

.200 

.128 

.512 

236 

.78o 

.100 

.400 

1.011 

.353 

1.120 

237 

.222 

.2<3 

.022 

.20C 

.128 

.524 

238 

.458 

.076 

.404 

1.749 

.529 

1.460 

239 

.22o 

.171 

.040 

.345 

.192 

.721 

24O 

.537 

.076 

.-73 

1.749 

.529 

1.53C 

241 

.307 

.171 

.054 

.345 

.192 

.73C 

242 

.391 

.130 

.204 

1.325 

.353 

.770 

.".43 

.442 

.130 

.298 

1.3-5 

.353 

.790 

244 

.4o3 

.130 

.312 

1.325 

.353 

.800 

245 

.592 

.100 

.682 

2.274 

.529 

1.344 

246 

1.024 

.223 

.234 

.449 

.192 

.635 

247 

.004 

.100 

.696 

2.274 

.529 

1.354 

248 

1.25e 

.223 

.233 

.449 

.192 

.665 

249 

.557 

.158 

.4o3 

1.617 

.353 

.731 

250 

.541 

.158 

.450 

1.017 

.353 

.723 

251 

.573 

.158 

.476 

1.617 

.353 

.738 

252 

.627 

.120 

.882 

2.756 

.529 

1.242 

253 

1.362 

.270 

.382 

.544 

.192 

.537 

254 

.oS7 

.120 

.904 

2.75o 

.529 

1.295 

255 

1.543 

.270 

.432 

.544 

.192 

.609 

2  56 

.211 

.056 

.337 

3.179 

.710 

1.430 

2  57 

.194 

.100 

.099 

1.011 

.353 

.794 

258 

.082 

.154 

.017 

.415 

.202 

.598 

2  59 

.ISO 

.048 

.,05 

4.320 

.885 

1.735 

26C 

.139 

.086 

.:97 

1.373 

.441 

.993 

2ol 

.070 

.134 

.020 

.563 

.253 

.724 

2o2 

.191 

.043 

.529 

5.502 

1.062 

2.120 

2o3 

.179 

.076 

.153 

1.749 

.529 

1.204 

2  04 

.001 

.118 

.022 

.718 

.303 

.850 

265 

.235 

.073 

.497 

4.170 

.710 

1.220 

266 

.05u 

.130 

.034 

1.325 

.353 

.627 

267 

.211 

.063 

.60  3 

5.630 

.885 

1.510 

263 

.113 

.112 

.104 

1.739 

.441 

.778 

269 

.Oi.5 

.lo8 

.01c 

.704 

.,53 

.5o5 

270 

.20 

.0  56 

.734 

7.153 

1.062 

1.800 

271 

.181 

.0*9 

.2  10 

2.274 

.529 

.960 

272 

.005 

.154 

.031 

.933 

.303 

.660 

273 

.231 

.089 

.601 

5.097 

.710 

1.070 

274 

.106 

.158 

.089 

1.617 

.353 

.544 

275 

.206 

.076 

.714 

6.793 

.e85 

1.340 

276 

.133 

.135 

.155 

2.159 

.441 

.675 

277 

.202 

.068 

.894 

e.671 

1.062 

1.663 

278 

.196 

.120 

.280 

2.756 

.529 

.843 

279 

.llo 

.187 

.009 

1.131 

.303 

.565 

285 

l.e55 

.288 

1.047 

1.697 

.045 

.500 

286 

1.786 

.257 

1.401 

1.513 

.045 

.508 

287 

1.568 

.207 

.978 

1.217 

.045 

.520 

288 

1.366 

.161 

.058 

.949 

.045 

.559 

289 

2.941 

.190 

.415 

.28C 

.C18 

.404 

290 

3.045 

.450 

1.960 

1.202 

.120 

.575 

291 

3.082 

.413 

1.486 

.913 

.109 

.526 

292 

3.143 

.392 

.910 

.554 

.087 

.425 

293 

3.093 

.300 

.684 

.432 

.087 

.482 

294 

3.151 

.191 

.429 

.270 

.087 

.654 

300 

1.919 

.401 

1.263 

1.189 

.009 

.380 

301 

1.814 

.395 

.996 

1.019 

.069 

.392 

302 

1.650 

.323 

.721 

.832 

.069 

.413 

303 

1.810 

.289 

.421 

.453 

.054 

.381 

304 

2.086 

.242 

.228 

.I08 

.036 

.392 

305 

2.523 

.400 

.322 

.247 

.168 

.413 

306 

3.972 

.3CC 

.225 

.113 

.131 

.510 

